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ABSTRACT 


A  long  standing  concern  of  the  Navy  has  been  the  need  to  reduce 
structural  vibrations  in  plates,  such  as  ship's  hulls.  Recently,  it  was 
proposed  to  use  waveguide  absorbers  as  a  means  to  reduce  structural 
vibrations.  A  waveguide  absorber  is  a  device  which  is  made  from  damping 
material  but  mounted  to  the  structure  at  a  point  rather  than  the  full 
surface.  A  waveguide  absorber  removes  and  then  dissipates  vibrational 
energy  from  the  structure  through  traveling  waves.  The  performance  of  a 
waveguide  absorber  is  determined  from  the  driving  point  impedance,  the 
ratio  of  force  to  velocity  at  the  attachment  point.  This  study  has  determined 
the  theoretical  driving  point  impedance  for  viscoelastic  circular  plates. 

The  constitutive  relations  for  viscoelastic  material  were  evaluated 
and  it  was  proven  that  a  complex  Young's  modulus  and  shear  modulus 
can  be  used  to  represent  the  viscoelastic  material.  Mindlin's  theory  for 
elastic  isotropic  plates,  with  complex  moduli,  was  used  to  solve  for  the 
driving  point  impedance.  The  solutions  are  in  terms  of  Bessel  functions. 
Poisson's  boundary  conditions  were  used  at  the  free  edge  of  the  circular 
plate,  and  clamped  boundary  conditions  were  used  at  the  attachment  point, 
which  is  at  the  center  of  the  plate.  A  non-numeric  language,  REDUCE, 
was  used  to  solve  for  the  driving  point  impedance  given  the  appropriate 
equations  of  motion  and  boundary  conditions.  The  complex  Bessel 
functions,  needed  for  the  results  of  the  REDUCE  program,  were  generated 
by  a  FORTRAN  program.  To  validate  the  results  of  this  study,  elastic  and 


viscoelastic  plates  were  tested  to  determine  their  driving  point  impedance. 
A  comparison  of  the  theoretical  results  and  the  experimental  results  shows 
that  there  is  agreement  for  all  cases  studied. 
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I.  INTRODUCTION 


A.  BACKGROUND 

The  reduction  of  structural  vibration  in  ships  and  submarines  has 
been  a  long  standing  concern  of  the  Navy.  Reducing  structural  vibrations 
can  increase  the  life  of  equipment  as  well  as  decrease  the  amount  radiated 
noise  into  the  water.  One  of  the  most  widely  used  methods  of  passive 
structural  vibration  reduction  has  been  2  addition  of  damping  material 
directly  to  the  surface  of  the  structure[Ref.  1  :  pp.  312  -362].  Recently  Ungar 
and  Kurzweil  [  Ref.  2  ]  have  proposed  the  use  of  a  waveguide  absorber  to 
reduce  structural  vibrations.  A  waveguide  absorber  is  a  damping  device 
that  is  attached  indirectly  to  the  structure  through  a  mount,  at  a  point, 
rather  than  directly  to  the  surface.  The  purpose  of  the  waveguide  absorber 
is  to  remove  energy  from  the  structure  and  then  dissipate  the  energy  as 
traveling  waves.  The  waveguide  absorber  is  expected  to  remove  more 
energy  per  unit  mass  of  damping  material  than  a  surface  damping 
treatment  can. 

Past  studies  of  waveguide  absorbers  have  been  by  Ungar  and 
Kurzweil  [Ref.  3] ,  by  Ungar  and  Williams  [Ref.4],  and  by  Lee  [Ref.  5]  . 
These  studies  show  that  a  waveguide  absorber  is  able  to  reduce  the 
magnitude  of  the  vibrations  in  a  structure.  The  theoretical  results  are 
based  upon  the  driving  point  impedance,  the  ratio  of  force  to  velocity  at  the 
attachment  point.  There  has  not  been  good  agreement  between  the 
experimental  results  and  the  theoretical  results  in  these  studies.  However, 
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in  the  study  by  Lee  [Ref.  5]  of  viscoelastic  beams,  it  was  determined  that 
using  an  underlining  theory  that  included  the  shear  deformation  and 
rotatory  inertia,  i.e.  the  Timoshenko  beam  theory,  produced  better  results 
than  those  theories  which  did  not  include  them,  i.e.,  the  Euler  beam  theory. 

B.  PURPOSE 

Based  upon  the  experimental  success  of  beam  waveguide  absorbers  to 
reduce  structural  vibrations,  there  is  reason  to  believe  that  a  two 
dimensional  geometry  may  provide  increased  damping  per  unit  mass  of 
damping  material  than  the  one  dimensional  beam  does.  The  purpose  of 
this  study  is  to  investigate  the  feasibility  of  using  viscoelastic  circular  plates 
as  waveguide  absorbers,  by  determining  the  driving  point  impedance  of  the 
plate.  The  driving  point  impedance  will  be  evaluated  over  a  frequency 
range  of  10  to  2000  Hz.  This  range  of  frequencies  was  chosen  because  data 
gathering  techniques  and  material  information  are  not  reliable  below  five  to 
ten  hertz,  and  the  distance  that  vibrations  or  noise  can  travel  in  water  above 
2000  Hz  is  very  small  and  normally  ignored.  It  is  the  goal  of  this  study  that 
the  driving  point  impedance  be  presented  an  easy  to  use  equation,  or  an 
easy  to  use  computer  program.  The  driving  point  impedance  should  be  able 
to  be  calculated  with  a  limited  number  of  inputs. 

C.  FORMAT 

This  study  has  been  broken  down  into  four  sections  even  though  the 
actions  occurred  at  the  same  time.  The  sections  are  theoretical  results, 
experimental  results,  a  comparison  of  the  theoretical  and  experimental 


results,  with  conclusions,  and  recommendations  for  future  studies.  The 
body  of  this  study  was  written  based  upon  all  of  the  plates  used,  but  only 
shows,  in  the  figures,  the  results  for  one  plate  size.  The  figures  for  the 
other  plate  sizes  are  found  in  Appendix  C.  This  was  done  so  as  to  keep  the 
body  of  this  study  as  small  as  possible. 


IL  THEORY 


A.  VISCOELASTIC  MATERIAL  MODELING 

In  the  past  study  by  Bailey  and  Chen  [Ref.  6],  it  was  shown  that  one  of 
the  best  methods  of  modeling  viscoelastic  material  is  by  separating  the 
elastic  and  viscoelastic  parts.  The  elastic  portion  is  not  time  dependent, 
where  as  the  viscoelastic  portion  is  time  dependent.  The  model  of 
viscoelastic  material  can  be  created  by  using  a  resistive  element  as  shown 
in  Figure  1  and  a  damping  element  as  shown  in  Figure  2,  where  a  is  the 
stress,  e  is  the  strain,  E  is  the  Young's  modulus,  and  C  is  the  damping  of 
the  material.  These  models  are  based  upon  electrical  impedance  where 
strain  is  the  equivalent  of  voltage,  stress  the  equivalent  of  current,  and  the 
inverse  of  the  characteristics  is  the  equivalent  of  impedance. 


Figure  L  Resistive  Element  ModeL 
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Figure  2.  Damping  Element  ModeL 


The  stress  versus  strain  relationship  for  the  resistance  elements  is: 

a 

e=E  (2.1) 


and  for  the  damping  element  is: 


(2.2) 


The  combination  of  a  resistive  element  and  a  damping  element  is  known  as 
the  Maxwell  model  [Ref.  7].  The  model  shown  in  Figure  3  is  based  on  a  two 
element  Maxwell  model,  where  the  full  element,  with  damping,  represents 
the  viscoelastic  portion  and  the  degenerate  element,  no  damping, 
represents  the  elastic  portion  of  the  material.  The  notation  used  is,  Ee  is 
the  elastic  Young's  modulus,  Ey  is  the  viscoelastic  Young's  modulus,  Cy  is 
the  viscoelastic  damping,  Ce  is  the  elastic  stress,  ay  is  the  viscoelastic 
stress,  a  is  the  total  stress,  and  e  is  the  total  strain  in  the  material. 
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For  the  viscoelastic  portion  the  strain  verse  stress  relationship  is: 


or  taking  the  time  derivative  of  equation  (2.5a)  gives: 


(2.5a) 


e  = 


(2.5b) 


It  is  concluded  that  the  total  stress  of  equation  (2.3)  can  be  rewritten  as: 


o  =  Eee+cv  (2.6) 

where  the  viscoelastic  stress  must  obey  the  relationship  of  equation  (2.5b). 
The  stress  versus  strain  relationship  can  be  characterized  as  the  normal 
elastic  relationship  plus  an  added  term  that  incorporates  the  viscoelastic 
effects. 

For  reference  purposes  the  coordinate  system  of  Figure  4  is  used  for 
the  discussion  that  follows. 


Figure  4.  Rectangular  Coordinate  System  for  a  Plate. 


For  thin  plates,  it  is  normal  to  assume  that  the  plate  undergoes  a  biaxial 
stress  state  [Ref.  8],  but  for  this  study  triaxial  stress  is  used  with  the 
modification  that  the  strain  in  the  z  direction  is  zero  and  that  the  stress  in 
the  z  direction  is  very  small  and  can  be  ignored.  This  assumption  will 
allow  the  inclusions  of  the  shear  deformation  and  rotatory  inertia  later. 

The  stresses  acting  a  small  piece  of  material  are  shown  in  Figure  5,  where 
oi,  or  oii,  is  the  normal  stress  on  the  i  face  in  the  i  direction,  x  jj  is  the  shear 
stress  on  the  i  face  in  the  j  direction,  which  is  equal  to  t  jj. 


x 


Figure  5.  Stresses  in  a  Cube  of  Material. 


Based  upon  the  model  in  Figure  3  with  resultant  equation  (2.6),  the  three 
dimensional  stress  versus  strain  relationships  for  isotropic  material 
become: 


ez  =  °  (2.7a) 

E 

a  x  =  — ^(e  x  +  v  eey  )  +  ox  y  (2.7b) 

1  ~  ve 

E 

°y=  1_*2(£y+Veex)+0yv  (2.7c) 

A  ve 

Ee 

x*y  “  2(1+  ve)Y*y  +  x*y  y  (2  7d) 

Ee 

X«~  2(1  + ve)Y«  +  X“v  (2.7e) 

Efi 

V  2(i  +  ve)Yy*  +  ty*v  (2.7f) 


where  ei  is  the  strain  in  the  i  direction,  y  jj  is  the  rotation  in  the  i-j  plane, 
which  is  equal  to  y jj,  and  the  sub-sub-script  of  v  is  for  the  viscoelastic 
portion.  The  viscoelastic  terms  must  also  meet  the  following  conditions: 


^  x  u  x  1 

A  V  A  V  X  /  • 

Ev  +  Cv  ~  1  _v2(Ex  +  VvEy) 

(2.8a) 

°yv  .  °yv  1  / 2.  ...  z.  \ 

E  +  C  2  (ey +  Vvex) 

v  1  —  v* 

(2.8b) 

x*y  v  x*yv  1 

Ev  +  Cv  “2(l  +  vv)Y*y 

(2.8c) 

X»v  Xxzv  1 

Ev  +  Cv  2(1  +  vv) Yx* 

(2.8d) 

xy*v  xy*v  1 

Ev  +  Cv  "2(l  +  vv)Yy* 

(2.8e) 
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The  use  of  these  constitutive  relations  proves  to  be  quite  cumbersome 
in  all  but  the  most  simplest  of  problems.  This  is  because  simultaneous 
differential  equations  must  be  solved  to  relate  the  stresses  to  strains.  By  the 
use  of  a  complex  model,  where  the  time  dependent  terms  are  replaced  by 
complex  terms,  some  of  these  problems  can  be  overcome.  Figure  6  shows 
the  Laplace  Transform  [Ref.  9],  of  the  model  shown  in  Figure  3. 


Figure  6.  Laplace  Transform  of  Viscoelastic  Material  ModeL 

By  applying  circuit  analysis  methods  to  the  model  shown  in  Figure  6,  an 
equivalent  model,  shown  in  Figure  7,  with  a  complex  Young's  modulus, 
E  ,  can  be  determined.  The  viscoelastic  elements  are  in  series  with  each 
other  and  together  are  in  parallel  with  the  elastic  element.  Solving  for  E* 


E*  Ee  (ev  +  8Cv) 

Rearranging  terms  and  substituting  jco  for  s,  where  j=  V-  1  and  co  is  the 
rotational  frequency  in  rads/sec,  in  equation  (2.9)  yields: 


(2.9) 


.  EE2+co2C2E  +co2C2E  coC  E2 

E*  ev  v  v  v  e  .  .•  v  v 

_ - - - - - +  J 


2  o«2 

Ev  +  0)TV 


2  o  2 

Ev  +  co2Cv 


The  resultant  complex  model  is  shown  in  Figure  7. 


e 


Figure  7.  Complex  Viscoelastic  Material  ModeL 


(2.10) 


The  complex  model  not  only  allows  the  use  all  of  the  elastic  theory,  with  its 
well  defined  mathematics,  but  has  the  added  benefit  that  most  of  the 
manufactures  of  viscoelastic  material  provide  material  data  in  this  form 
[Ref.  1,  pp.  363-440].  Instead  of  leaving  the  complex  Young's  modulus  in 
terms  of  the  elastic  and  viscoelastic  parameters,  as  in  equation  (2.10),  three 
new  terms  are  defined  and  used.  The  real  part,  E',  is  defined  as: 


,  EE2+o)2C2E  +g)2C2E 

t/  _  e  v  v  v  v  e 

**  2  o  2 

Ev  +  ©2Cv 


(2.11) 
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The  imaginary  part,  E",  by: 


Ef  = 


(oC  E2 

V  V 
2  9 

Ev  +  £D2C 


2 

v 


The  ratio  of  the  imaginary  part  to  the  real  part  is  known  as  the  loss 
damping  factor,  rj,  and  defined  as: 


(2.12) 


The  complex  Young's  modulus  may  now  be  written  as: 


(2.13) 


E*  =  Ef  +  jEf 


(2.14) 


or  as: 


E  =Ef(l+jTi)  (2.15) 

It  is  normal  to  assume  that  the  Poisson's  ratio.v  ,  for  the  viscoelastic 
material  is  the  same  for  the  elastic  and  viscoelastic  portions,  which  leads  to 
the  result  that  the  shear  modulus  is  also  a  complex  parameter,  G*,  and  is 
related  to  the  complex  Young's  modulus  by: 


G*  = 


E 


2  (1  +  v) 


=  G  +jG 


(2.16) 


The  real  and  imaginary  parts,  respectively,  of  the  Young's  modulus  and 
shear  modulus  are  related  in  identical  manner.  The  loss  factor  is  the  same 
for  both  moduli. 


12 


Elastic  theory  can  now  be  used  to  analyze  structures  using 
viscoelastic  material  when  Young's  modulus  and  shear  modulus  are 
replaced  by  the  viscoelastic's  complex  moduli.  This  is  a  great  benefit  in  the 
analysis  of  structures  with  viscoelastic  plates  because  a  separate  and  new 
governing  theory  does  not  have  to  be  created. 

B.  ELASTIC  SOLUTION 

For  this  study  an  elastic  thin  plate  theory  that  indudes  the  effects  of 
shear  deformation  and  rotatory  inertia  was  needed  to  determine  the  driving 
point  impedance  of  the  viscoelastic  plate.  The  use  of  such  a  theory,  for 
beams,  in  the  study  by  Lee  [Ref.  5]  provided  better  agreement  between 
theoretical  and  experimental  results  than  did  theories  which  did  not 
incorporate  these  effects.  Mindlin  [Ref.  10]  has  presented  such  a  theory  for 
isotropic  thin  plates.  This  work  showed  that  isotropic  elastic  thin  plates 
can  be  described  by  a  set  of  three  simultaneous  partial  differential 
equations.  It  is  based  upon  assuming  that  the  displacements  in  the  x 
direction,  u,  and  in  the  y  direction,  v,  of  Figure  4,  vary  linearly  in  the  z 
direction  and  the  displacement  in  the  z  direction,  w,  is  not  a  function  of  z, 
or  written  in  equation  form: 


u  =  z-  yx(x,  y,t) 

(2.17a) 

v=z-  yy(x,  y,  t) 

(2.17b) 

w  =  W(x,  y,  t) 

(2.17c) 
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where  \|/x  anf^  Vy  are  potential  functions  to  be  determined.  This  implies 
that  the  displacements  in  the  x  and  y  directions,  u  and  v,  vary  as  shown  in 
Figure  8. 


Figure  8.  Variations  in  z  and  y  displacements  as  a  function  of  z 
for  the  Minrflin  Plate  Umary. 

Mindlin's  original  work  was  in  a  rectangular  coordinate  system,  but 
was  adapted  to  cylindrical  coordinate  system  by  Mindlin  and 
Deresiewicz  [Ref.  11]  and  then  updated  by  Irie,  Yamada,  and 
Aomura  [Ref.12],  and  to  curvilinear  by  Callahan  [Ref.13].  The  notation  of 
Irie,  Yamada,  and  Aomura  [Ref  12]  is  used  in  this  study  owing  to  studying 
circular  plates  and  their  use  of  dimensionless  parameters.  The  reference 
directions  of  the  cylindrical  coordinate  system  are  shown  in  Figure  9. 
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r 


Figure  9.  Cylindrical  Coordinate  System  for  a  Plate. 

The  displacements,  in  terms  of  cylindrical  coordinates,  which  are 
equivalent  to  those  of  equation  (2.17)  are  the  radial  displacement,  ur,  the 
angular  displacement,  uq,  and  the  displacement  in  the  z  direction,  w. 

They  are  expressed  mathematically  as: 


ur  =  z-  yrt(r,  9,  t) 


(2.18a) 


ua  =  z  •  yet(r,  6,  t) 


(2.18b) 


w  =  Wt(r,  0,  t) 


(2.18c) 


where  Vrt  and  V0t  are  potential  functions  that  are  to  be  determined.  The 

moments  and  shears  in  a  plate  are  shown  in  Figure  10  and  defined  in 
equation  (2.19)  [Ref.  10]  for  the  case  of  cylindrical  coordinate  system. 
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where  Mii,  or  Mi,  is  the  bending  moment  on  the  i  face  per  unit  of  length  in 
the  r-6  plane,  Mjj,  which  is  equal  to  Mji,  is  the  twisting  moment  on  the  i 
face  per  unit  of  length  in  the  r-0  plane,  and  Qi,  is  the  shear  on  the  i  face  per 
unit  of  length  in  the  r-0  plane  .  The  plate  geometry  is  shown  in  Figure  8 
and  Figure  9. 


Figure  UL  Plan  and  Side  View  of  Plate  Showing  Geometry. 

For  the  case  of  an  isotropic  plate  in  cylindrical  coordinates,  the 
Mindlin  governing  equations  are  [Ref.ll]: 


17 


(2.20a) 


jMg  m„-mm  ph3  aVr 

dr  r  do  r  12  ^2 

, i**»  , 2M^  Q  _ph!^ 
dr  +r  90  r  ^0  12  ^2 


9Qr  1 aQe 
Or  +  r  ae 


(2.20b) 


(2.20c) 


where  p  is  the  density  of  the  material.  It  is  not  desirable  to  solve 
equations  (2.20)  in  the  present  form,  because  of  the  simultaneous 
differential  equations,  but  rather  in  a  simpler  form.  Mindlin  [Ref.  10] 
showed  that  the  three  simultaneous  second  order  partial  differential 
equations,  can  be  combined  to  produce  a  single  fourth  order  partial 
differential  equation,  in  terms  of  w.  The  only  drawback  of  using  this  single 
equation  is  that  it  can  not  handle  the  full  number  of  geometric  boundary 
conditions.  Another  suggestion  by  Mindlin  [Ref.  10]  was  to  assume  that  the 
displacements  were  all  periodic  functions  of  time  of  the  type  expfcot).  By 
doing  this,  the  three  simultaneous  equations  can  be  converted  into  three 
independent  second  order  partial  differential  equations,  which  are  much 
easier  to  solve,  but  still  use  of  all  the  geometric  boundary  conditions. 

For  this  study,  it  was  desired  to  use  the  geometric  boundary 
conditions  directly,  in  order  to  improve  the  agreement  between  the 
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theoretical  and  experimental  results,  and  since  most  shipboard  structural 
vibrations  are  from  rotating  equipment,  periodic  in  nature,  the  final 
suggestion  by  Mindlin  will  be  used.  As  such,  the  displacements  are 
assumed  to  be  periodic  of  the  form  expfcot ),  and  expressed  as: 


Ur  =  z-  vr(r,  0)ej(Ot 

(2.21a) 

u0  =  z-  Ve(r,  0)ej<Bt 

(2.21b) 

w=W(r,0)eja)t 

(2.21c) 

where  the  unknowns  are  yr,  yg,  and  W.  It  is  convenient  to  define  the 

following  dimensionless  parameters  that  are  based  upon  the  plate  geometry 

and  the  frequency  [Ref.  12]: 

R_(i)2 

R-  12 

(2.22a) 

s  D  1  Ch'>2 

i^Ga^  6(1 

(2.22b) 

^4  pha4to2 
*  "  D 

(2.22c) 
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(2.22d) 


52=K 

l  2 


R+S  +  {(K-S)2  +  -i-l 

1  x4J 


2  1,4 

52"  2l 


“ 

1' 

2 

R  +  S-< 

(r-s>2  +  4| 

. 

l  X  J  J 

(2.22e) 


52  = 


3  1  -  V 


(2.220 


c  = 


1 


(2.22g) 


5 


c  = 


2  4  i 

^-1 


(2.22h) 


where  D  is  the  flexural  rigidity  expressed  by: 


D  = 


Eh' 


12(l  -v2) 


(2.23) 


and  k  is  the  shear  coefficient  taken  from  Mindlin  [Ref.10].  The  shear 
coefficient  is  used  to  take  into  account  that  the  shear  is  not  a  constant  value 
as  assumed,  but  varies  across  the  thickness  of  the  plate.  For  elastic  plates 
the  shear  coefficient  is  normally  set  equal  to: 


k2=^ 

12 


(2.24) 


For  the  presently  assumed  case  of  periodic  motion  of  isotropic  plates, 
the  governing  partial  differential  equations  become  [Ref.  12]: 


(v2+<k=° 

(2.25a) 

(v2+822)w2=0 

(2.25b) 

(V2+S3>3  =  ° 

(2.25c) 

where  V2  is  the  Laplacian  operator  and  the  variable  wi,  w2,  and  W3  are 
dummy  variables  which  are  related  to  the  unknowns,  \yr,  Yq,  and  W  by  the 

following  relationships  [Ref.  12]: 


dw  dw  ,  dw 

v  r  “  (<Y - x) sr + (°  2 - x)~sr  +  rle- 


jdWj 


1*”2 


ve  =  (°1  -  J)  r le"  +  (°2  '  J)  rle  5T 


W  =  Wj  +  w2 


(2.26a) 


(2.26b) 


(2.26c) 
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It  is  much  easier  to  work  with  the  governing  equations  as  shown  in 
Equation  (2.25)  than  those  of  Equation  (2.20).  The  problem  has  been  reduced 
from  having  to  solve  simultaneous  partial  differential  equations,  to  solving 
three  independent  partial  differential  equations.  Equation  (2.25)  is  a  form  of 
Bessel's  partial  differential  equation  which  has  a  solution  in  terms  of 
Bessel  functions  [Ref.  14]  [Ref.15,  pp.  22-26].  For  this  study  the  solution  to 
equation  (2.25)  is  assumed  to  be: 

Wl  =  [V"(8ia)  +A4Yn(81i)]cos(ne)  (2.27a) 

W2  =  [A2Jn(S2?)  +A5Yn(S2i)>°8(D0)  (2  27b) 

W3  =  [A3*^n(83a)  +  A  <227c) 

where  Ai  are  constants  to  be  set  by  the  boundary  conditions,  Jn (x)  is  the 
Bessel  function  of  the  first  kind  of  order  n,  and  Yn(x)  is  the  Bessel  function 
of  the  second  kind  of  order  n.  The  Bessel  function  of  the  second  kind  of 
order  n,  has  been  added  to  the  solution  of  Irie,  Yamada,  and  Aomura 
[Ref. 12]  because  of  the  hole  in  the  center  of  the  plates  being  used  in  this 
study. 

Six  boundary  conditions  are  needed  to  completely  establish  the 
solution  to  this  plate  problem.  For  the  plates  under  study,  three  boundary 
conditions  are  at  the  center  of  the  plate,  at  r=b,  and  three  boundary 
conditions  are  at  the  edge  of  the  plate,  at  r=a.  The  boundary  conditions  are 
given  in  terms  of  displacements  and  the  stress  in  the  plate.  At  the  center, 


r=b,  for  a  plate  being  used  as  a  waveguide  absorber,  it  must  be  mounted  to  a 
fixture.  It  has  been  assumed  that  the  fixture  is  rigid,  has  a  radius  of  b,  and 
clamps  the  plate  such  that  there  is  no  twisting  and  that  the  displacement  in 
the  z  direction  at  r=b  is  the  same  for  the  fixture  as  the  plate.  In 
mathematical  terms,  the  boundary  conditions  at  the  plate's  center  are: 


Vr(b,0)=O 

(2.28a) 

¥0(b.  0)  =0 

(2.28b) 

W(b,6)=w0 

(2.28c) 

If  the  plate  was  not  attached  to  the  fixture,  or  if  b=0,  then  the  constants  A4, 
A5,  and  Ag  would  have  to  be  set  to  zero  to  insure  that  the  solution  remained 
bounded,  as  the  Bessel  functions  of  the  second  kind  are  unbounded  as  the 
argument  approaches  zero. 

At  the  free  edge  of  the  plate,  r=a,  the  boundary  conditions  are  no 
bending  moments,  twisting  moments,  or  shearing  on  the  free  edge,  or  in 
mathematical  terms: 

Mrr(a,6)  =  0  (2.29a) 

Mrt(a,6)  =  0 


(2.29b) 


These  are  known  as  Poisson's  boundary  conditions  [Ref.16].  Poisson's 
boundary  conditions  can  be  used  because  the  governing  differential 
equations  were  not  reduced  to  form  a  single  differential  equation  which 
ignores,  or  reduces  the  effects  of  the  shear  deformation  and  rotary  inertia. 
Had  a  reduced  differential  equation  been  used,  then  Kirchoff  s 
boundary  [Ref.  16]  [Ref.17]  conditions  would  have  to  be  used.  Kirchoff  s 
boundary  conditions  for  the  free  edge  are: 


Mr(a,0)=O 


(2.30a) 


Qr(r, 


e)  + 


1  dM^r,  0)- 


do 


=  0 


(2.30b) 


r  =  a 


Kirchoff  s  boundary  conditions  are  an  attempt  to  include  the  three 
geometric  boundary  conditions  of  Poisson  that  exists  at  the  free  edge  into 
two  equations.  The  necessity  for  this  arose  from  a  plate  theory,  equivalent 
to  the  Euler  beam  theory,  which  could  not  incorporate  three  edge  boundary 
conditions  but  only  two.  The  present  theory,  Mindlin’s,  does  not  give  up 
boundary  conditions  for  simplicity.  The  use  of  a  reduced  theory  and 
Kirchoff  s  boundary  condition  can  lead  to  a  large  disparity  between  the 
theoretical  and  the  experimental  displacements  [Ref.  18]  when  r 
approaches  the  free  edge  at  a. 
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All  that  is  left  to  allow  analyzing  the  plate  is  to  relate  the 
displacements,  or  potential  functions,  to  the  moments  and  shears.  The 
equations  below  accomplish  this  [Ref.  12] : 


(2.31a) 

Mee  =  d{f(v,  +  ^r)+ 

(2.31b) 

(2.31c) 

Qr  =  K»Gh(yr  +  ^) 

(2.31  d) 

(2.31  e) 

Finding  the  constants  Ai  to  Aq,  is  a  straight  forward  matter  of 
combining  equation  (2.27)  into  equation  (2.26),  and  then  into  equation  (2.31) 
from  which  the  boundary  conditions  of  equations  (2.28)  and  (2.29)  can  be 
satisfied.  Once  the  constants  Ai  to  Ag  are  known,  the  potential 
functions, Yr  and  yq,  the  displacements,^,  ug,  and  w,  and  the  moments 
and  shears,  Mrr,  Mqq,  Mrg,  Qr»  and  Qg,  are  known  anywhere  in  the  plate. 

With  this  information  the  driving  point  impedance  can  be  calculated. 


Mechanical  impedance,  Z,  is  defined  as: 


Z 


(2.32) 


where  F  is  the  force  at  the  point  of  interest  and  V  is  the  velocity  at  the  point 
of  interest.  For  this  study,  the  point  of  interest  is  the  the  driving  point 
where  the  fixture  mounts  to  the  plate.  Based  upon  the  assumption  made 
about  the  fixture  and  nature  of  the  displacements,  the  force,  F,  at  the 
driving  point  can  be  calculated  by: 

2n 

f=  JrQ^b.eje^ae  (23. 

0 


and  the  velocity,  V,  is  related  to  the  displacement  at  the  driving  point  by: 


tt  a(w(b,0)eja>t)  .  „T/1L  jtot  .  j«* 

V  =  — - - ^ - =jcoW(b,  0)e  =j©wQe 


(2.34) 


Combining  equations  (2.33)  and  (2.34)  into  equation  (2.32)  yields: 


2k 

JrQr(b,e)ae 

_0 _  (2.35) 

jG>w0 


Equation  (2.35)  is  the  driving  point  impedance. 


The  steps  necessary  to  solve  for  the  driving  point  impedance  have 
now  been  laid  out.  The  only  difficulty  lies  in  determining  the  constants  Ai 
to  A6  for  any  plate  geometry  (a,b,h)  and  material  (p,E*,v).  The  work 
necessary  to  determine  these  constants  is  beyond  the  capability  of  normal 
hand  written  calculations.  Assuming  that  the  driving  point  impedance 
can  be  derived,  the  driving  point  impedance  must  be  calculated  over  and 
over  for  different  frequencies,  material,  and  geometry.  The  only  way  to 
solve  for  the  driving  point  impedance  is  with  the  aid  of  a  computer.  The  use 
of  the  computer  was  split  into  two  parts.  The  first  was  to  use  a  symbolic 
language  to  algebraically  determine  the  driving  point  impedance,  and  the 
second  was  to  use  a  numerical  language  to  calculate  the  driving  point 
impedance  for  a  given  set  of  conditions. 

C.  COMPUTER  IMPLEMENTATION 

Figure  12  shows  the  flow  path  used  in  the  symbolic  program  to 
determine  the  driving  point  impedance.  One  of  the  primary  considerations 
in  the  selection  of  a  non  numerical  language,  which  will  be  used  to 
determine  the  driving  point  impedance,  was  the  ease  of  differentiating 
Bessel  functions.  The  language  of  REDUCE  [Ref.  19}  has  the  capability  to  be 
easily  taught  how  differentiate  functions,  so  was  used  in  this  study. 
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Methods  for  Differentiation 
of  Bessel  Functions 


Figure  12.  Symbolic  Language  Program  Flow  Chart. 


The  REDUCE  program  did  the  differentiation  of  the  Bessel  functions 
and  substitutions  as  required  very  well.  The  program  showed  that  the 
resultant  motion  was  symmetric  with  n  equal  to  zero.  This  agrees  with  the 
boundary  conditions  imposed  at  the  center  of  the  plate,  motion  only  in  the  z 
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direction  with  no  twisting.  It  was  only  at  the  point  of  trying  to  determine 
the  constants  Ai  to  A6  did  REDUCE  have  trouble.  In  effect,  the  REDUCE 
program  was  solving  a  linear  six  by  six  symbolic  matrix.  In  the  present 
implementation  of  REDUCE  on  the  Naval  Postgraduate  School’s  IBM  3033 
computer  there  was  insufficient  storage  to  to  find  all  six  constants  directly. 
In  order  to  alleviate  this  problem,  the  program  was  broken  up  into  two 
parts,  one  part  to  set  up  the  necessary  equations  to  be  solved,  and  the  other 
to  solve  the  simultaneous  linear  equations.  Cramer's  rule  of  linear  matrix 
solution  [Ref.  20]  was  employed  to  solve  the  simultaneous  linear  equations 
through  the  use  of  dummy  variables.  Even  using  this  approach  there  was 
insufficient  memory  to  solve  the  system  of  equations.  The  symmetry  of  the 
problem  was  then  taken  advantage  of,  A3  and  Aq  do  not  need  to  be  solved 
because  sin(nO)  is  zero.  The  result  was,  that  the  REDUCE  program  only 
had  to  solve  a  four  by  four  matrix  for  Ai ,  A2,  A4,  and  A5,  which  it  was  able 
to  do.  Once  these  constants  were  found  then  the  magnitude  of  the  force,  F, 
magnitude  of  the  velocity,  V,  and  the  driving  point  impedance,  Z,  were 
found.  The  driving  point  impedance  is  in  terms  of  Bessel  functions 
evaluated  at  the  center,  r=b,  and  free  edge,  r=a,  of  the  circular  plate.  As  a 
final  step  in  the  REDUCE  program  the  driving  point  impedance  was 
written  in  FORTRAN  type  statements  to  allow  a  simpler  generation  of  a 
FORTRAN  program.  The  REDUCE  programs  used  in  this  study  are  in 
Appendix  A. 

The  flow  chart  for  the  numerical  program  used  to  calculate  the 
driving  point  impedance  is  shown  in  Figure  13.  The  numerical  language  of 
FORTRAN  was  used  in  this  study. 
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The  FORTRAN  program  was  written  to  so  as  to  minimize  the 
number  of  inputs  required  in  order  to  calculate  the  driving  point  impedance 
over  a  range  of  frequencies.  It  was  decided  that  the  inputs  would  be  the  type 
of  material,  geometry,  and  temperature  of  the  plate,  and  the  upper  and 
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lower  frequency  bounds  of  interest.  The  remaining  information  needed  to 
calculate  the  driving  point  impedance  would  either  be  calculated  or  stored 
in  table  faction. 

As  was  noted  by  Irie,  Yamada,  and  Aomura  [Ref.  12],  the  parameters 
62^  and  832  can  become  negative  such  that  the  argument  of  the  Bessel 

function  is  a  complex  value.  In  their  study,  to  get  around  this,  whenever 
82 2  or  832  become  negative,  the  assumed  solution  went  from  the  Bessel 

function  of  the  first  kind  to  the  modified  Bessel  function  of  the  first  kind  and 
the  argument  of  the  Bessel  function  was  made  into  a  real  quantity.  When 
dealing  with  viscoelastic  plates,  where  the  argument  of  the  Bessel  functions 
is  always  complex,  this  strategy  was  not  acceptable.  It  was  pointed  out  by 
McLachlan  [Ref.  15,  p.  27]  that  the  Bessel  functions  of  the  first  and  second 
kind  are  still  solutions  to  the  Bessel's  partial  differential  equation  even  with 
complex  arguments.  This  means  that  the  solution  can  be  calculated  if  the 
Bessel  functions  are  in  complex  notation.  All  of  the  Bessel  function 
algorithms  found  on  the  Naval  Postgraduate  School's  IBM  3033  are  written 
for  real  valued  arguments,  not  complex  arguments,  so  could  not  be  used.  A 
computer  algorithm  had  to  be  written  to  find  the  values  of  the  Bessel 
function  of  the  first  kind,  second  kind,  and  their  first  and  second 
derivatives.  It  was  decided  that  the  simplest  method,  and  most  straight 
forward,  though  not  the  most  elegant,  would  be  to  write  the  algorithm  with 
the  Taylor  series  expansion  of  the  functions.  It  must  also  be  pointed  out 
that  this  algorithm  was  written  for  the  special  case  of  Bessel  function  order, 
n,  being  zero.  The  Taylor  series  expansions  used  to  calculate  the  Bessel 
functions  and  their  first  and  second  derivatives  are  [Ref  21]: 
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(2.36b) 


(2.36c) 


(2.36d) 


(2.36e) 


(2.360 


where  x  is  a  dummy  variable  and  y  is  Euler's  number  which  is  equal  to: 


T‘„^-(1  +  2+S+-  •+-5-‘n<n>) 


(2.37) 
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For  this  study  Eider's  number  was  approximated  by: 


y*  0. 57721 5664901 5328606065:  (2.38) 

The  code  used  to  calculate  the  Bessel  functions  was  written  in  quadruple 
precision  complex  FORTRAN,  which  is  an  extension  of  FORTRAN  77  found 
on  the  Naval  Postgraduate  School's  IBM  3033  computer.  Using  this 
precision  of  code  for  the  Taylor  series  expansion  insured  that  the  resultant 
values  had  an  accuracy  of  at  least  six  significant  places. 

In  order  to  validate  the  programs  written  for  this  study,  the  driving 
point  impedance  of  elastic  plates  would  also  be  calculated  and  compared 
experimentally.  This  did  not  present  any  problems  in  the  use  of  the 
programs,  because  an  elastic  material's  Young'  modulus  and  shear 
modulus  can  be  represented  by  complex  values  with  zero  imaginary  part. 
The  only  addition  required  to  the  programs  is  the  elastic  material's 
characteristics.  For  the  elastic  material  aluminum  was  chosen,  and  for 
the  viscoelastic  material  LD-400  was  chosen. 

It  has  been  pointed  out  by  Mindlin  [Ref.10]  ,and  Reismann  and 
Pawlik  [Ref.  22]  that  the  value  of  shear  coefficient,  k  in  equation  (2.24),  is  not 
fixed  but  depends  on  the  material's  Poisson's  ratio  and  the  frequency  of  the 
first  antisymmetric  mode  of  thickness-shear  vibration.  Equation  (2.24)  is 
based  upon  only  the  first  antisymmetric  mode  of  thickness-shear  vibration 
and  corresponds  to  a  Poisson's  ratio  of  v  =  0.176.  For  both  the  elastic 
material  and  viscoelastic  material  under  study  the  Poisson's  ratio  is 
greater  than  this.  It  was  decided  that  two  values  of  the  shear  coefficient 
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will  be  used  for  each  type  of  material  being  studied.  The  values  used  are 
shown  in  Table  1 . 

The  material  Young's  modulus,  or  shear  modulus,  and  the 
Poisson's  ratio  must  also  be  supplied.  For  the  aluminum  plate  this 
information  was  taken  from  Timoshenko  [Ref.  23].  For  the  viscoelastic 
plates,  the  information  was  taken  from  Lee  [Ref.  5] ,  which  was  based  upon 
data  from  Nashif,  Jones,  and  Henderson  [Ref.  1,  pp.  390*392],  and  from 
material  data  sheet  supplied  by  the  manufacture,  United  McGill  [Ref.  24]. 
Only  the  data  from  United  McGill  [Ref.  24]  was  temperature  compensated 
in  the  program  written  for  this  study.  The  information  from  United 
McGill  [Ref.  24]  was  provided  on  a  log-log  scale  chart  which  allowed 
determination  of  complex  Young’s  modulus  at  three  different 
temperatures,  68,  86,  and  104  degrees  Fahrenheit.  For  each  temperature 
the  complex  Young's  modulus  was  determined  at  four  frequencies,  10, 100, 
1000,  and  10,000  Hz.  The  information  was  taken  in  terms  of  the  log  of  the 
value,  so  that  the  data  was  better  conditioned  for  use  in  the  numerical 
calculations.  A  four  point,  third  order  polynomial  curve  fit  [Ref.  25]  was 
applied  at  each  of  the  above  temperatures  to  allow  calculating  the  modulus 
at  other  frequencies.  At  plate  temperatures  between  those  taken,  a  linear 
interpolation  was  used.  The  polynomial  equations  used  are  shown  below. 
For  a  temperature  of  68  degrees  Fahrenheit: 
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D  „  =  2. 025  -  0. 06250 f,  +  0.  lOOOf?  -  0.  01250f?  (2  39a) 

E  log  log  log  v  ' 

D  =  1.775 -0.07084 f  +  0.1125f?  -0.01667ff  (2.39b) 

E"  log  log  log 

For  a  temperature  of  86  degrees  Fahrenheit: 

D  =1.000  +  0. 6125f.  -  0. 1875ff  +0.02500ff  (2.40a) 

E  log  log  log 

D  E„  =  1 . 050  +  0.  3875flog  -  0. 1  OOOf^  +  0.  01250^  (2.40b) 

For  a  temperature  of  104  degrees  Fahrenheit: 

D  =0.9000  +0.  08750f  +0.01250ff  +  0.  0000001 590  ff  (2.41a) 

E  log  log  log 

D  =  0. 8000  +  0.  2542 f  -  0. 06250ff  +  0.  008333 ff  (2.41b) 

E  log  log  log 


where  flog  is  related  to  the  frequency  by: 

f,  =  log  (f) 
log 


(2.42) 


DE'  and  De"  are  dummy  variables  which  are  related  to  the  real  part,  E', 
and  the  imaginary  part,  E",  of  the  complex  modulus  by: 
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(2.43a) 


E'  =  10 


E"  =  10 


(de) 

(de) 


1450.5 
1450. 5 


(2.43b) 


where  E’  and  E"  are  in  the  units  of  pounds  per  square  inch.  The  material 
characteristics  for  the  three  reference  temperatures  are  shown  in  Figures 
14  to  16.  The  linear  temperature  interpolation  was  done  in  terms  of  De'  and 
De"  and  then  converted  into  E’  and  E"  by  equation  (2.43).  It  was  felt  that 
aluminum  would  be  insensitive  to  temperature  changes  and  the  data  from 
Nashif,  Jones,  and  Henderson  [Ref.  1,  pp.  390-392]  was  too  coarse  to  allow 
accurate  calculations  other  than  at  room  temperature,  75  degrees 
Fahrenheit.  The  shear  modulus  and  loss  factor  equations  are  repeated 
here  for  the  data  from  Nashif,  Jones,  and  Henderson  [Ref.  1,  pp.  390-392]: 

G'  =  0. 00002503  f  3  -  0. 1752  f  2  +  457. 5883  f  +  29280  (2.44a) 


Tl  = 


-  Q  5273 


0. 65e 1 


(2.44b) 


The  moduli  and  loss  factor  for  the  viscoelastic  material,  from  both 
sources  are  shown  in  Figures  17  to  19  for  a  temperature  of  75  degrees 
Fahrenheit.  This  shows  for  the  same  material,  there  is  little  agreement  on 
the  material  characteristics  by  different  sources. 
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Figure  14.  Viscoelastic  Material  Characteristics  from 

United  McGill  [Ref.  24]  at  68  Degrees  Fahrenheit 
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Figure  15.  Viscoelastic  Material  Characteristics  from 

United  McGill  [Ret  24]  at  86  Degrees  Fahrenheit. 
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Figure  18.  Comparison  of  Viscoelastic  Material  Shear  Modulus  from 
Nashifj  Jones,  Henderson  [Ret  1,  pp.  390-392]  and 
United  McGill  [Ret  24] . 
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Figure  19.  Comparison  of  Viscoelastic  Material  Loss  Factor  from 
Nashif,  Jones,  Henderson  [Ref.  i,  pp.  390*392]  and 
United  McGill  [Ref  24] . 
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Other  inputs  to  the  FORTRAN  program  are  the  density  and  the 
temperature  of  the  plate.  The  summary  of  the  constants  used  are  listed  in 
Table  1 .  A  copy  of  the  FORTRAN  programs  used  can  be  found  in 
Appendix  B.  The  driving  point  impedance  was  calculated  at  100  point 
across  the  frequency  range  of  interest.  For  the  case  of  the  elastic  plate,  the 
theoretical  impedance  at  a  mode  will  tend  to  infinity.  However,  for  the 
frequency  range  of  interest,  100  to  2000  Hz,  the  chances  of  having  the 
frequency  of  interest  being  the  same  frequency  as  a  mode  are  very  small, 
and  an  overflow  of  the  computer  is  unlikely.  The  results,  for  all  materials, 
are  shown  in  the  real  and  imaginary  parts  of  the  driving  point  impedance 
due  to  the  waveguide  predictions  being  based  upon  the  real  part  of  the 
impedance.  In  an  effort  to  check  the  calibration  of  the  experimental  portion 
of  this  study,  the  driving  point  impedance  was  calculated  in  the  range  of  10 
to  105  Hz.  The  lower  limit,  10  Hz,  was  set  by  the  material  characteristics  of 
United  McGill  [Ref.  24]  not  going  less  than  this,  while  the  upper  limit 
corresponded  to  the  lower  limit  of  the  normal  frequency  range  with  a  5  Hz 
overlap. 


TABLE  L  DIMENSIONS  AND  PROPERTIES  OF 
MATERIALS  UNDER  STUDY. 


ELASTIC 

VISCOELASTIC 

MATERIAL 

AL 

LD-  400  - 

-DM-375 

Reference 

Ref.  23 

Ref.  1 

Ref.  24 

h  (in) 

.2- 

8 

a  (in) 

4,  5,  and  6 

b  (in) 

JL 

8 

E  (psi) 

11.0  M 

E'  (psi) 

Equations 

(2.39)  to 
(2.43) 

E"  (psi) 

G'  (psi) 

Equation 

(2.44) 

n 

V 

0.33 

0.5 

*2 

_ 

K2 

12 

and 

0.8600 

_  2 

JL-  and  0.9128 

X  « 

3.06713  „ 
jtlO-3 

1.71101  _ 
xlO-3 

1.70736  „ 
*10-3 
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The  results  for  the  elastic  plates,  shown  in  Figures  20  and  21 ,  show 
that  the  plate  is  acting  as  a  lumped  mass  in  the  low  frequency  range  of  10  to 
105  Hz.  A  lumped  mass  has  no  real  part  to  the  impedance,  and  the 
imaginary  increases  at  a  linear  rate,  which  is  what  is  happening  here. 

The  results  for  the  high  frequency  range,  100  to  2000  Hz,  do  not  provide 
much  information,  except  the  frequency  of  the  first  mode,  as  shown  in 
Figures  22  and  23. 

The  results  for  the  viscoelastic  plates,  shown  in  Figures  24  to  27, 
show  that  there  is  a  significant  difference  in  the  driving  point  impedance 
due  to  the  material  characteristics  used,  either  Nashif,  Jones,  and 
Henderson  [Ref.  1,  pp.  390-392]  or  United  McGill  [Ref.  24].  It  also  shows 
that  at  very  low  frequencies  the  viscoelastic  material  is  behaving  as  a 
lumped  mass,  with  little  real  part  and  a  linearly  increasing  imaginary 
part.  The  effect  of  the  temperature  compensation  of  the  material 
characteristics  supplied  by  United  McGill  [Ref.  24]  is  demonstrated  by  the 
results  shown  in  Figures  28  and  29.  This  shows  that  the  impedance 
decreases  and  becomes  smother  with  increases  in  temperature. 

In  all  cases  the  effect  of  a  different  shear  coefficient,  k,  is 

« 

insignificant  and  a  shear  coefficient  of  2  (k=0.9069)  will  be  used  for 

the  remainder  of  this  study. 
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The  dominant  factors  for  the  viscoelastic  plates  which  effect  the 
calculated  driving  point  impedance,  are  the  source  of  material 
characteristics  and  plate  temperature  for  a  fixed  plate  geometry,  and  the 
geometry  of  the  plate.  Figures  74  to  77  and  Figures  96  to  99,  in  Appendix  C, 
show  the  results  for  the  five  and  four  inch  plates. 
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Effect  of  Shear  Coefficient  on  the  Theoretical  Real  Part  of 
the  Driving  Point  Impedance  for  a  6  in  Radius  Elastic 
Plate  in  the  Frequency  Range  of  10  to  106  Hz. 
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Figure  2L 
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Effect  of  Shear  Coefficient  on  the  Theoretical  Imaginary 
Part  of  the  Driving  Point  Impedance  for  a  6  in  Radius 
Elastic  Plate  in  the  Frequency  Range  of  10  to  105  Hz. 


(Nl/baS-itTl)  33MVQHdKI 


Figure  22.  Effect  of  Shear  Coefficient  on  the  Theoretical  Real  Part  of 
the  Driving  Point  Impedance  for  a  6  in  Radius  Elastic 
Hate  in  the  Frequency  Range  of  100  to  2000  Hz. 
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Figure  23.  Effect  of  Shear  Coefficient  on  the  Theoretical  Imaginary 
Part  of  the  Driving  Point  Impedance  for  a  6  in  Radius 
Elastic  Plate  in  the  Frequency  Range  of  100  to  2000  Hz. 


(Ni/oas-ier)  aoMvgaawi 


Figure  24.  Effect  of  Shear  Coefficient  and  Material  Characteristics  on 
the  Theoretical  Beal  Part  of  the  Driving  Point  Impedance 
for  a  €  in  Radius  Viscoelastic  Plate  in  the  Frequency 
Range  of  10  to  105  Hz  at  a  Temperature  of  75.0  Deg.  F. 
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Figure  25.  Effect  of  Shear  Coefficient  and  Material  Characteristics  on 
the  Theoretical  Imaginary  Part  of  the  Driving  Point 
Impedance  for  a  6  in  Radius  Viscoelastic  Plate  in  the 
Frequency  Range  of  10  to  105  Hz  at  a  Temperature  of  75.0 
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Figure  26.  Effect  of  Shear  Coefficient  and  Material  Characteristics  on 
the  Theoretical  Real  Part  of  the  Driving  Point  Impedance 
for  a  6  in  Radius  Viscoelastic  Plate  in  the  Frequency 
Range  of  100  to  2000  Hz  at  a  Temperature  of  745  Deg:  F. 
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Figure  28.  Effect  of  Shear  Coefficient,  Material  Characteristics,  «nH 
Temperature  on  the  Theoretical  Real  Part  of  the 
Driving  Point  Impedance  for  a  6  in  Radius  Viscoelastic 
Plate  in  the  Frequency  Range  of  100  +o  2000  Hz. 


in.  EXPERIMENT 


A.  EXPERIMENT  SET  UP 

A  total  of  6  plates  were  tested  during  this  study.  Three  were  made 
out  of  aluminum,  an  elastic  material,  and  three  were  made  out  of  LD-400 
material,  a  viscoelastic  material.  All  were  3/8  of  an  inch  thick  with  a  plate 
from  each  material  at  a  radius  of  four,  five,  and  six  inches  as  listed  in 
Table  1 .  The  viscoelastic  plates  where  made  from  twelve  inch  square  tiles 
that  are  3/8  inch  thick,  this  limited  the  maximum  plate  size  to  a  six  inch 
radius.  The  other  two  plates  were  chosen  so  that  there  would  be  a  definite 
change  in  characteristics  for  a  change  in  geometry  and  so  that  the  ratio  of 
thickness  to  radius  did  not  become  too  large  and  deviate  from  thin  plate 
theory.  To  determine  the  driving  point  impedance  of  the  plates,  the 
equipment  and  configuration  of  Figure  30  was  used. 


In  Figure  30,  WR  is  for  Wilcoxon  Research  and  HP  for  Hewlett  Packard. 

The  Wilcoxon  Research  F4/F7  shaker  and  Z7  measurement  head  are  an 
integral  unit  that  was  mounted  on  a  stand  for  the  impedance 
measurements.  The  driving  signal  to  the  shaker  is  provided  by  the  Hewlett 
Packard  HP-3562A  dynamic  signal  analyzer  via  the  Wilcoxon  Research 
PA7  power  amplifier  and  N7/N9  matching  network.  The  force  and 
acceleration  measurements  of  the  plate  were  provided  by  the  Wilcoxon 
Research  Z7  impedance  head  mounting  base,  and  then  amplified  by  a 
Kistler  503  charge  amplifier,  for  the  force  signal,  or  by  a  PCB  462 A  charge 
amplifier,  for  the  acceleration  signal,  before  going  to  the  HP-3562A  dynamic 
signal  analyzer.  The  HP-3562A  dynamic  signal  analyzer  was  operated  in 
the  swept  sine  mode  for  all  of  the  tests  to  insure  that  the  coherence  of  the 
data  was  at,  or  almost  at,  one.  A  minimum  of  400  sample  points  were 
taken  when  the  frequency  range  of  interest  was  100  to  2000  Hz,  and  a 
minimum  of  160  sample  points  were  taken  when  the  frequency  range  of 
interest  was  5  to  105  Hz.  These  were  only  a  minimum  as  the  analyzer  was 
set  to  provide  auto  frequency  resolution,  or  take  more  sample  points,  when 
the  incoming  data  showed  a  large  change  from  the  previous  data  point. 

The  analyzer  was  also  set  up  to  control  the  amount  of  power  going  to  the 
shaker  so  that  the  incoming  signals  were  of  a  quality  that  would  keep  the 
coherence  at  one.  In  the  frequency  range  of  less  than  2000  Hz,  the  data  was 
taken  with  the  F7  piezoelectric  vibration  generator  electrically  secured  and 
all  data  gathered  with  the  F4  electromagnetic  vibration  generator.  This 
was  done  because  the  F7  was  producing  considerable  ring  in  the  force  and 
acceleration  signals,  which  was  degenerating  the  quality  of  the  signal 
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going  to  the  analyzer.  The  F7  is  normally  used  for  the  high  frequency 
measurements  and  the  F4  for  the  low  frequency  measurements  [Ref.26].  A 
photograph  of  the  four  inch  radius  plate  mounted  on  the  shaker  is  shown  in 
Figure  31 . 
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The  easiest  method  of  determining  the  driving  point  impedance  of  the 


plate  from  the  information,  force  and  acceleration,  being  provided  to  the 
analyzer,  is  to  start  with  the  plate's  frequency  response  and  then 
mathematically  convert  this  information  into  the  driving  point  impedance 
by  math  functions  built  into  the  analyzer.  The  frequency  response,  H(f) ,  as 
measured  by  the  analyzer,  is: 


H(f)  = 


GFA(f> 

GFF^G 


(3.1) 


where  GFA(f)  i8  the  cross  spectrum  between  the  force  signal,  F,  and  the 
acceleration  signal,  A,  and  GFF(0  is  the  power  spectrum  of  the  force  signal 
[Ref.  27].  The  transfer  function,  H(f),  is  the  ratio  of  acceleration  to  force  or: 

H<f>  -  f  (3.2) 

The  acceleration  is  related  to  the  velocity  by: 


A  = 


(3.3) 


By  combining  equations  (3.3)  and  (3.2)  and  rearranging  the  resultant 
equation  the  driving  point  impedance  becomes: 


jca 

"  H(f) 


(3.4) 


In  the  analyzer  the  frequency  response,  H(f),  was  first  multiplied  by  (jto)*l 
and  then  the  reciprocal  of  the  resultant  function  was  taken  to  be  the  driving 
point  impedance,  Z. 

B.  PLATE  MOUNTING  DESIGN 

Some  way  of  attaching  the  plate  to  the  shaker  must  be  designed  and 
used.  For  use  as  a  waveguide  absorber,  the  plate  will  have  a  hole  drilled  at 
the  center  of  the  plate  for  a  bolt  or  screw.  This  idea  was  used  as  the  starting 
point  for  the  design  of  the  plate  to  shaker  mount.  A  mount  was 
manufactured  from  steel  rod  two  inches  long  with  3/8  inch  threads  on  one 
end,  the  shaker  end,  and  1/4  inch  threads  on  the  other  end,  the  plate  end. 
The  plate  was  held  between  two  1/4  inch  nuts  and  washers  and  a  lock 
washer  as  shown  in  Figure  32. 


63 


1 

4 


J  3"  1. 
1  8  1 


Figure  32.  Steel  Plate  to  Shaker  Mount 


This  mount  proved  to  be  unsatisfactory  in  operation.  Its  weight  and  size 
were  totally  against  it.  It  resonated  right  in  the  middle,  from  about  700  to 
1200  Hz,  of  the  frequency  range  of  interest,  100  to  2000  Hz.  The  data  taken, 
when  using  this  mount,  was  useless.  A  new  mounting  system  was 
designed,  using  aluminum,  to  minimize  the  weight  and  size  jf  the  mount 
and  is  shown  in  Figure  33. 
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Attach  to  Shaker 


Figure  33.  Aluminum  Plate  to  Shaker  Mount 

This  mount  was  extremely  light  weight  and  has  almost  no  effect  on  the 
plate  data.  There  was  sufficient  friction  produced  in  the  attachment  to  the 
shaker,  that  the  plate  and  mount  did  not  vibrate  lose  during  the  data  runs, 
and  a  lock  washer  was  not  required.  As  a  means  of  checking  the 
effectiveness  of  this  mount,  the  driving  point  impedance  of  the  shaker  and 
mount,  with  no  plate,  were  determined  over  the  same  frequency  range  as 
with  the  plate.  For  these  runs,  the  analyzer  had  the  same  configuration  as 
with  plates  except  that  the  auto  frequency  resolution  was  turned  off.  If  the 
auto  frequency  resolution  were  left  on,  the  data  runs  would  have  taken  too 
long,  greater  than  2  hours.  The  results  of  these  runs  showed  that  the 
shaker  and  mount  were  acting  as  if  a  pure  mass  of  1.48656  mslugs.  The 
shaker  and  mount  mass  is  much  less  than  the  mass  of  the  smallest  plate 
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under  study,  so  its  effects  are  ignored.  A  3/8  inch  hole  was  drilled  in  the 
center  of  each  plate  to  allow  mounting.  The  plate  sat  on  a  5/8  inch  outside 
diameter  boss  on  the  shaker,  and  then  was  clamped  down  by  the  aluminum 
mount.  Both  the  shaker  and  mount  have  an  outside  diameter  of  5/8  inch 
which  was  used  in  all  of  the  theoretical  calculations  for  the  value  of  inner 
radius,  b. 

C.  VERIFICATION  OF  IMPEDANCE  SCALE 

In  all  of  the  data  runs,  the  calibrated  output  from  the  force  and 
acceleration  transducer,  plus  the  gain  of  the  charge  amplifiers,  was 
inputted  into  the  analyzer  to  insure  that  the  results  of  the  analyzer  were  of 
the  proper  scale.  As  a  means  of  checking  the  scale  of  these  inputs,  the 
driving  point  impedance  of  an  elastic  plate  can  be  compared  experimentally 
and  theoretically  at  a  certain  frequency.  At  frequencies  less  than  the  first 
mode,  the  elastic  plate  can  be  considered  as  a  lumped  mass  as  shown  in 
Figure  34. 


Mass,  m 


Figure  34  Lumped  Mass  Model  of  Elastic  Plate. 
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Slimming  the  forces  on  the  mass,  results  in  the  equation  of  motion: 


F  -  ma  =  F  -  mx  =  0 


(3.5) 


The  acceleration,  a,  in  terms  of  velocity,  V,  is: 

9V  . 

a  -  -5T  =J“V 


(3.6) 


The  force,  F,  in  terms  of  the  velocity,  V,  is: 


F  =  jiomV 


(3.7) 


The  impedance,  Z,  then  becomes: 

Z  =  y  =jcom 


(3.8) 


This  allows  calculating  the  impedance  of  the  plate,  based  upon  its  weight, 
or  mass,  and  the  frequency  of  interest.  A  frequency  of  100  Hz  was  chosen 
for  this  comparison  on  the  elastic  plates  under  study.  The  100  Hz  frequency 
is  less  than  the  lowest  mode  of  the  plates  under  study.  The  results  of  this 
comparison  as  shown  in  Table  2.  Only  the  imaginary  parts  of  the 
impedance  are  looked  at,  as  the  model  shows  there  is  no  real  part,  only  an 
imaginary  part. 
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TABLE  2.  IMAGINARY  PART  OF  IMPEDANCE  FOR 
ELASTIC  PLATES. 


Data  Run 

Plate 

Radius 

Calculated 

5-105  Hz 

100-2K  Hz 

Impedance,  Imaginary  Part,  lbf-sec/in 

4  in 

3.1178 

2.77836 

2.79759 

5  in 

4.8973 

4.22329 

4.45218 

6  in 

7.0535 

6.79937 

6.9406 

The  results  of  this  comparison  show  that  the  analyzer  and  the  associated 
equipment  have  the  proper  impedance  scale. 

D.  DATA  RUNS 

Each  of  the  six  plates  used  in  this  study,  were  subjected  to  two 
different  frequency  ranges  at  room  temperature.  One  data  run  was  from  5 
to  105  Hz,  for  the  purpose  of  checking  the  operation  and  calibration  of  the 
equipment,  and  the  other  from  100  to  2000  Hz,  the  frequency  range  of 
interest.  As  an  added  measure  both  of  the  six  inch,  elastic  and  viscoelastic, 
plates  were  run  from  100  to  2000  Hz  at  an  elevated  temperature  of  about  100 
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degrees  Fahrenheit  to  check  the  effect  of  temperature  on  Young'  modulus. 
The  surface  temperature  of  the  plate  was  measured  with  an  Omega  871 
digital  thermometer  with  a  type  K  surface  temperature  probe.  A  number  of 
points  on  the  plate  were  measured  before  and  after  the  data  run.  It  was 
observed  that  the  shaker  was  adding  heat  to  the  plates  during  the  data  run. 
This  resulted  in  a  temperature  gradient  of  about  five  degrees  Fahrenheit 
from  the  center  of  the  plate  to  the  outer  edge.  The  temperature  of  the  plate 
was  taken  as  the  average  of  all  of  the  readings. 

On  all  data  runs  the  coherence  remained  at  one  for  every  sample 
point  taken.  The  repeatability  of  the  data  was  tested  on  six  inch  radius 
elastic  and  viscoelastic  plates  by  repeating  the  data  runs  after  being 
removed  from  the  shaker  and  then  reinstalled,  and  by  running  the  swept 
sine  run  in  the  decreasing  mode  instead  of  the  increasing  mode.  There 
was  almost  no  difference  in  the  driving  point  impedance  for  these  runs  and 
the  normal  data  run.  It  is  concluded  that  the  data  is  repeatable. 

E.  EFFECT  OF  THE  SHAKER  ON  THE  IMPEDANCE 
MEASUREMENTS 

The  shaker  and  the  mount  have  mass  which  contributes  to  the  data 
measurements.  The  effect  of  this  mass  can  be  taken  into  account  by 
considering  the  model  of  the  measurement  system  shown  in  Figure  35. 
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Figure  35.  Impedance  Measurement  ModeL 


The  measured  impedance  for  this  model  is: 

Z  ■  Z 

7  =7  ,  —Mount - Plate 

Measured  Shaker  Z  +  Z  '*■ 

Mount  Plate 

If  the  impedance  of  the  mount  can  be  ignored,  then  the  plate  impedance 
becomes  the  measured  impedance  minus  the  shaker  impedance  or: 


Z  —  z  —  z 

Plate  Measured  Shaker 


(3.10) 
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This  calculation  can  be  accomplished  in  the  analyzer  very  easily.  In  order 
to  determine  the  impedance  of  the  shaker  and  the  mount,  data  runs  were 
taken  of  the  shaker  and  mount  together,  and  with  the  shaker  alone.  The 
results  of  these  data  runs,  Figures  36  to  39,  show  that  the  shaker  and 
mount  are  effectively  a  lumped  mass  with  little  real  impedance  and  a 
linearly  increasing  imaginary  impedance  with  frequency.  This  occurred 
over  both  frequency  ranges,  5  to  105  Hz  and  100  to  2000  Hz,  tested.  The  mass 
of  the  mount  was  measured  to  be  0.34946  mslugs,  while  the  mass  of  the 
shaker,  from  m=Z/co,  is  1.0374  mslugs.  The  mass  of  the  lightest  plate  is 
30.0123  mslugs,  which  is  almost  100  times  greater  than  the  mount  mass 
and  30  times  greater  than  the  shaker  mass.  As  such  the  mass  effect  of  the 
mount  has  been  ignored  and  the  plate  impedance  found  by  equation  (3.10). 

It  was  noted  that  the  data  below  about  15  Hz  was  very  hard  to  accurately  get 
without  the  mass  of  a  plate  being  present.  The  data  in  this  region  has  a  fair 
amount  of  scatter. 
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Figure  37. 


Experimental  Imaginary  Part  of  the  Impedance  of  Shaker, 
With  and  Without  Mount,  in  the  Frequency  Range  of  5  to 
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Figure  38.  Experimental  Real  Part  of  the  Impedance  of  Shaker,  With 
and  Without  Mount,  in  the  Frequency  Range  of  100  to 
2000  Hz. 
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Figure  39.  Experimental  Imaginary  Part  of  the  Impedance  of  Shaker, 
With  and  Without  Mount,  in  the  Frequency  Range  of  100 
to  2000  Hz. 
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In  Figures  40  to  43  the  difference  between  the  measured  and 
corrected  plate  impedance  is  shown.  There  is  some  small  difference  in  the 
case  of  imaginary  impedance  and  no  difference  for  the  real  impedance,  the 
shaker  appears  as  a  lumped  mass.  However  the  difference  between  the 
measured  and  corrected  impedance  is  not  considered  to  be  significant  to 
the  conclusions  of  this  study,  and  as  such,  the  data  will  be  in  the 
uncorrected  form  for  the  rest  of  this  study. 
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Figure  40.  Experimental  Real  Part  of  the  Driving  Point  Impedance  of 
a  6  in  Radius  Viscoelastic  Plate,  With  and  Without  Shaker 
and  Mount  Correction,  in  the  Frequency  Range  of  5  to  105 
Hz  at  a  Temperature  of  75.0  Deg.  F. 


Figure  4L  Experimental  Imaginary  Part  of  the  Driving  Point 

Impedance  of  a  6  in  Radius  Viscoelastic  Plate,  With  and 
Without  Shaker  and  Mount  Connection,  in  the  Frequency 
Range  of  5  to  106  Hz  at  a  Temperature  of  75.0  Deg.  F. 
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Figure  42.  Experimental  Real  Part  of  the  Driving  Point  TmppAtnfw  of 
a  6  in  Radius  Viscoelastic  Plate,  With  and  Without  Shaker 
and  Mount  Correction,  in  the  Frequency  Range  of  100  to 
2000  Hz  at  a  Temperature  of  74J5  Deg.  F. 
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Figure  43.  Experimental  Imaginary  Part  of  the  Driving  Print 

Impedance  of  a  6  in  Radius  Viscoelastic  Plate,  With  and 
Without  Shaker  and  Mount  Correction,  in  the  Frequency 
Range  of  100  to  2000  Hz  at  a  Temperature  of  74.5  Deg.  F. 
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F.  RESULTS  OF  DATA  RUNS 


For  the  elastic  plates  in  the  low  frequency  range  of  5  to  105  Hz,  shown 
in  Figures  44  and  45,  the  plates  appear  to  lumped  mass  as  expected,  with  a 
small  real  part  and  a  linearly  increasing  imaginary  part  with  frequency. 

In  the  normal  frequency  range,  100  to  2000  Hz,  the  results,  Figures  46  and 
47,  for  the  elastic  plates  hold  little  information,  only  showing  the  location  of 
the  first  mode,  as  was  the  case  in  the  theoretical  portion.  There  was  almost 
no  difference  in  the  driving  point  impedance  for  the  six  inch  elastic  plate 
when  data  was  taken  at  room  temperature  and  when  taken  at  an  elevated 
temperature,  as  shown  in  Figures  48  and  49.  This  confirms  that  Young's 
and  shear  moduli  for  the  aluminum  plate  need  not  be  temperature 
compensated. 


Figure  45.  Experimental  Imaginary  Part  of  the  Driving  Point 

Impedance  of  a  6  in  Radius  Elastic  Plate  in  the  Frequency 
Range  of  5  to  105  Hz. 
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Figure  46.  Experimental  Real  Part  of  the  Driving  Point  Impedance  of 
a  6  in  Radius  Elastic  Plate  in  the  Frequency  Range  of  100 
to  2000  Hz. 
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Figure  47.  Experimental  Imaginary  Part  of  the  Driving  Point 

Impedance  of  a  6  in  Radius  Elastic  Plate  in  the  Frequency 
Range  of  100  to  2000  Hz. 
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Figure  48.  Effect  of  Temperature  on  the  Experimental  Real  Part  of 
the  Driving  Point  Impedance  for  a  6  in  Radius  Elastic 
Plate  in  the  Frequency  Range  of  100  to  2000  Hz. 
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Figure  49.  Effect  of  Temperature  on  the  Experimental  Imaginary  Part 
of  the  Driving  Point  Impedance  for  a  6  in  Radius  Elastic 
Plate  in  die  Frequency  Range  of  100  to  2000  Hz. 


The  driving  point  impedance  for  the  viscoelastic  plate  in  the  low 
frequency  range  of  5  to  105  Hz  is  shown  in  Figures  50  and  51.  The  results 
for  the  normal  frequency  rage  of  100  to  2000  Hz  are  shown  in  Figures  52  and 
53.  There  is  a  major  effect  on  the  driving  point  impedance  by  the 
temperature  of  the  plate,  as  was  shown  in  the  theoretical  calculations. 

This  is  shown  in  Figures  54  and  55. 

Although  the  accuracy  of  the  data  is  effected  by  such  things  as  the 
calibration  of  equipment  and  the  actual  plate  geometry,  it  is  dominated  by 
the  material  characteristics  used.  The  material  characteristics  are  a 
design  specification  provided  by  the  vendor  and  the  material  may  have  a 
significant  statistical  variation  in  the  material  properties.  The  plate 
temperature  also  effects  the  material  conditions  for  a  fixed  plate  geometry 
and  a  different  temperature  measurement  technique  may  provide  a  slight 
increase  in  the  accuracy. 
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50.  Experimental  Real  Part  of  the  Driving  Point  Impedance 
for  a  6  in  Radius  Viscoelastic  Plate  in  the  Frequency 
Range  of  10  to  106  Hz  at  a  Temperature  of  75.0  Deg.  F. 


Figure  5L  Experimental  Imaginary  Part  of  the  Driving  Point 

Impedance  for  a  6  in  Radius  Viscoelastic  Plate  in  die 
Frequency  Range  of  10  to  106  Hz  at  a  Temperature  of  75.0 
Deg.  F. 
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Figure  62.  Experimental  Real  Part  of  the  Driving  Point  Impedance 
Cor  a  6  in  Radius  Viscoelastic  Plate  in  the  Frequency 
Range  of  100  to  2000  Hz  at  a  Temperature  of  745  Deg:  F. 
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Figure  53.  Experimental  Imaginary  Part  of  the  Driving  Point 

Impedance  for  a  6  in  Radius  Viscoelastic  Plate  in  the 
Frequency  Range  of  100  to  2000  Hz  at  a  Temperature  of 
74.5  Deg:  F. 


92 


(Nl/OaS-iai)  33MV(I3dM 

Figure  54.  Effect  of  Temperature  on  the  Experimental  Real  Fart  of 
the  Driving  Point  Impedance  for  a  6  in  Radios 
Viscoelastic  Plate  in  the  Frequency  Range  of  100  to 


IV.  COMPARISON 


A.  ELASTIC  PLATES 

For  the  elastic  plates  in  the  frequency  range  of  10  to  105  Hz  the 
agreement  between  the  theoretical  results  and  the  experimental  results  is 
very  good  as  shown  by  Figures  56  and  57.  In  the  frequency  range  of  100  to 
2000  Hz  the  only  thing  that  can  be  compared  is  the  location,  in  frequency,  of 
the  impedance  peak,  as  shown  in  Figures  58  and  59.  For  the  three  elastic 
plates  studied  the  experimental  peak  occurred  at  a  lower  frequency,  by 
about  100  Hz,  than  the  theoretical  peak.  In  the  high  frequency  range,  the 
theoretical  calculations  entirely  missed  the  spike  in  the  real  portion  of  the 
driving  point  impedance,  due  to  the  frequency  resolution  used  in  the 
calculations. 

It  is  concluded  that  the  ideas  and  programs  used  in  this  study  have 
been  validated  to  be  correct  by  these  results. 
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Figure  56.  Theoretical  and  Experimental  Comparison  of  the  Real  Part 
of  the  Driving  Point  Impedance  for  a  6  in  Radius  Elastic 
Plate  in  the  Frequency  Range  of  5  to  105  Hz. 
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Figure  57.  Theoretical  and  Experimental  Comparison  of  the 

Imaginary  Part  of  the  Driving  Point  Impedance  for  a  6  in 
Radius  Elastic  Plate  in  the  Frequency  Range  of  5  to 
105  Hz. 
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Figure  58.  Theoretical  and  Experimental  Comparison  of  the  Real  Part 
of  the  Driving  Point  Impedance  for  a  6  in  Radius  Elastic 
Plate  in  the  Frequency  Range  of  100  to  2000  Hz. 
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Figure  59.  Theoretical  and  Experimental  Comparison  of  the 

Imaginary  Part  of  the  Driving  Point  Impedance  for  a  6  in 
Radius  Elastic  Plate  in  the  Frequency  Range  of  100  to 
2000  Hz. 


99 


B.  VISCOELASTIC  PLATES 

In  all  cases  there  was  almost  no  agreement  between  the 
experimental  results  and  the  theoretical  results  when  the  material 
characteristics  of  Nashif,  Jones,  and  Henderson  [Ref.  1,  pp.  390-392]  were 
used.  But,  as  shown  by  Figures  60  to  63,  there  is  agreement  between  the 
experimental  results  and  the  theoretical  results  when  the  material 
characteristics  of  United  McGill  [Ref.  24]  are  used. 

The  remaining  discussion  deals  with  theoretical  results  with 
material  characteristics  supplied  by  United  McGill  [Ref. 24],  In  the  low 
frequency  range  the  experimental  results  would  break  or  change  at  a  lower 
frequency  than  the  theoretical  results  would.  This  effect  is  similar  to  the 
experimental  first  mode  for  the  elastic  plates  being  lower  than  the 
theoretical.  In  the  high  frequency  range  the  experimental  real  part  of  the 
driving  point  impedance  was  always  lower  than  the  theoretical,  while  the 
experimental  imaginary  part  of  the  driving  point  impedance  was  always 
higher  than  the  theoretical.  It  is  as  if  the  theoretical  calculations  are  a  few 
degrees  in  phase  off  from  the  experimental  results.  Figures  64  and  65  show 
the  importance  of  including  the  temperature  in  the  determination  of  the 
material  characteristics.  The  theoretical  and  experimental  results  have 
the  same  shape  but  are  shifted  by  a  small  amount.  It  is  felt  that  this  shift  is 
due  to  the  uncertainty  in  the  material  characteristics  used  in  the 
theoretical  calculations. 
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Theoretical  Experimental  Comparison  of  the  Real  Part 
of  the  Driving  Point  Impedance  for  a  6  in  Radius 
Viscoelastic  Plate  in  the  Frequency  Range  of  5  to  106  Hz  at 
a  Temperature  of  75.0  Deg.  F. 
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Theoretical  and  Experimental  Comparison  of  the  Real  Part 
of  the  Driving  Point  Impedance  for  a  6  in  Radios 
Viscoelastic  Hate  in  the  Frequency  Range  of  100  to  2000  Hz 
at  a  Temperature  of  745  Deg.  F. 
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Theoretical  and  Experimental  Comparison  of  the 
Imaginary  Part  of  the  Driving  Point  Impedance  for  a  6  in 
Radios  Viscoelastic  Plate  in  the  Frequency  Range  of  100  to 
2000  Hz  at  a  Temperature  of  745  Deg.  F. 
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Figure  64.  Theoretical  and  Experimental  Comparison  of  the  Real  Part 
of  the  Driving  Point  Impedance  for  a  6  in  Radius 
Viscoelastic  Hate  in  the  Frequency  Range  of  100  to  2000  Hz 
at  a  Temperature  of  100.5  Deg.  F. 
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Figure  66.  Theoretical  and  Kipwrimflnfail  Comparison  of  the 

Imaginary  Part  of  the  Driving  Point  Impedance  for  a  6  in 
Radius  Viscoelastic  Plate  in  the  Frequency  Range  of  100  to 
2000  Hz  at  a  Temperature  of  100i5  Deg.  F. 


C.  CONCLUSIONS 


1.  The  driving  point  impedance  of  a  circular  viscoelastic  plate 
can  be  predicted  by  the  programs  developed  by  this  study  by  knowing  the 
plate's  geometry  and  the  material  characteristics.  The  shape  of  the 
theoretical  and  experimental  results  is  the  same,  but  are  shifted  a  small 
amount  due  to  the  material  conditions  used. 

2.  The  effect  of  differing  shear  coefficients  on  the  driving  point 
impedance  is  immaterial  and  can  be  ignored. 

3.  The  driving  point  impedance  of  a  viscoelastic  plate  is  controlled 
by  the  plate  geometry  and  material  conditions.  The  material  conditions  are 
greatly  effected  by  the  temperature  of  the  plate. 

4.  For  viscoelastic  plates,  the  magnitude  of  the  real  and 
imaginary  part  of  the  driving  point  impedance  at  a  mode  are  about  the 


V.  RECOMMENDATIONS 


Based  upon  the  results  of  this  study  the  following  recommendations 
are  made  for  studies  continuing  in  this  area. 

1 .  Investigate  the  effect  of  twisting  of  the  waveguide  absorber  mount  on 
the  driving  point  impedance. 

2.  In  the  past  study  by  Lee  [Ref.  5]  a  constrained  layer  beam  waveguide 
absorber  was  used  with  some  success.  A  constrained  layer  circular 
plate  can  be  used  as  a  waveguide  absorber  and  its  characteristics  need 
to  be  studied. 

3.  The  effectiveness  of  the  present  viscoelastic  circular  waveguide 

absorbers  on  vibration  control  need  to  be  studied  and  compared  to  a 
surface  damping  treatment  and  a  point  mass,  all  of  the  same  weight. 
Investigate  whether  the  circular  plate  waveguide  absorbers  can 
provide  more  energy  dissipation  per  unit  mass  or  volume  of  damping 
device  than  other  forms  of  vibration  control? 

4.  This  study  has  shown  that  plate  geometry  and  plate  temperature 
effect  the  driving  point  impedance.  A  sensitivity  study  should  be 
undertaken  to  further  investigate  this  dependance  of  the  driving  point 
impedance  to  changes  in  plate  geometry  and  plate  temperature. 


5.  Investigate  the  ability  to  design  a  circular  waveguide  absorber  to 
remove  a  set  amount  of  energy,  at  a  point  on  the  structure,  at  a  certain 
frequency  band.  The  design  should  be  based  upon  the  driving  point 
impedance  of  the  waveguide  absorber  and  driving  point  impedance  of 
the  structure  at  the  attachment  point. 

6.  Investigate  the  ability  to  design  and  use  multiple  waveguide 
absorbers  to  dampen  over  a  wide  frequency  range  or  temperature 


range. 


APPENDIX  A. 

A.  REDUCE  PROGRAM  TO  SOLVE  A  FOUR  BY  FOUR  MATRIX. 


%  THIS  PROGRAM  IS  DESIGNED  TO  SOLVE  A  4  BY  4  MATRIX; 
K  MATRIX  METHOD; 

X  PROGRAM  TCA4A; 

OFF  NAT$ 

OUT  TT1H$ 

LINELENGTH  (72); 

D1=*DET  MATC(B11,B22,B33,B44),(C21,C22,C23,C24), 
(C31,C32,C33,C34), (C41 ,C42, C43, C44) ; 

D2 : =DET  MATCCC11,C12,C13,C14),CB11,B22,B33,B44), 
CC31,C32,C33,C34),(C41,C42,C43,C44); 

D3:=DET  MAT((C11,C12,C13,C14), CC21 , C22, C23, C24) , 
(B11,B22,B33,B44), ( C41 , C4? , C43 , C44 ) ; 

D4 : =DET  MAT((C11,C12.C13,C14), ( C21 , C22 , C23 , C24 ) , 
(C31,C32,C33,C34),  CB11,B22,B-33,B44); 

Di=DET  MAT ((C11,C12,C13,C14),(C21,C22,C23,C24), 

(C31,  C32,C33,C34) , (C41 ,C42, C43,C44) ; 

SHUT  TT1H$ 

BYE; 
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B.  REDUCE  PROGRAM  TO  DETERMINE  THE  DRIVING  POINT 


IMPEDANCE. 


XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX; 

%  PROGRAM  TA6C; 

X  SPECIAL  CASE  OF  N=0; 

OFF  NAT$ 

OUT  TT1F$ 

X  THIS  PROGRAM  ATTEMPTS  TO  FI HD  THE  SIX(6)  FACTORS  REQUIRED  TO; 

X  DESCRIBE  THE  WAVE  MOTION  OF  A  PLATE  WITH  FREE  EDGES  AT  R=A,  AND; 

X  A  DRIVEN  CLAMPED  CONDITION  AT  THE  CENTER  R=B . ; 

X  SPECIAL  CASE  OF  N=0; 

X  LINE  SIZE  TO  ALLOW  FILE  STORAGE; 

LINELEHGTH  (72); 

X  ADDED  DIFFERENTIATION  RULES  FOR  THE  BESSEL  FUNCTIONS; 

X  JN  OF  DELTAI ; 

OPERATOR  JN1, JN1DF1, JH1DF2, JN1DF3; 

FOR  ALL  R  LET  DF(JN1(R) , R ) =DELTA1XJN1DF1 ( R)/A; 

FOR  ALL  R  LET  DFC JN1 DF1 (R) , R) =DELTA1XJN1 DF2( R )/A ; 

FOR  ALL  R  LET  DFC JN1DF2C R) , R) =DELTAlx JN1DF3C R)/A; 

X  YN  OF  DELTAI; 

OPERATOR  YN1 , YN1DF1 , YN1DF2, YN1DF3; 

FOR  ALL  R  LET  DFC YN1 (R ) , R) =DELTA1XYN1DFI ( R)/A; 

FOR  ALL  R  LET  DFC YN1DF1 ( R) , R) =DELTA1XYN1DF2C R)/A; 

FOR  ALL  R  LET  DFCYN1DF2C R) , R) =DELTA1XYN1DF3C R)/A; 

X  JN  OF  DELTA2; 

OPERATOR  JN2, JN2DF1, JN2DF2, JN2DF3; 

FOR  ALL  R  LET  DFC JN2CR) ,R)=DELTA2XJN2DF1(R)/A; 

FOR  ALL  R  LET  DFC JN2DF1 ( R) , R) =DELTA2XJN2DF2( R)/A; 

FOR  ALL  R  LET  DFC JN2DF2C R) , R) *DELTA2*JN2DF3CR)/A; 

X  YN  OF  DELTA2; 

OPERATOR  YH2 , YN2DF1 . YN2DF2 , YN2DF3 ; 

FOR  ALL  R  LET  DFCYN2CR) , R)=DELTA2XYN2DF1CR)/A; 

FOR  ALL  R  LET  DFC YN2DF1 C R) , R) =DELTA2XYN2DF2C R)/A; 

FOR  ALL  R  LET  DFC YN2DF2C R) , R ) =DELTA2*YN2DF3( R)/A; 

X  JN  OF  DELTA3; 

OPERATOR  JN3 , JN3DF1 , JN3DF2 , JN3DF3; 

FOR  ALL  R  LET  DFC JN3C R ) , R ) =DELTA3*JN3DF1 C R)/A ; 

FOR  ALL  R  LET  DFC JN3DF1 ( R ) , R ) =DELTA3XJN3DF2( R )/ A; 

FOR  ALL  R  LET  DFC JN3DF2CR) ,R)=DELTA3*JN3DF3CR)/A; 

X  YN  OF  DELTA3; 

OPERATOR  YN3. YN3DF1 , YN3DF2, YN3DF3; 

FOR  ALL  R  LET  DFC YN3C R ) , R ) =DELTA3*YN3DF1 ( R)/ A ; 

FOR  ALL  R  LET  DF(YN3DF1(R),R)=DELTA3*YN3DF2(R)/A; 

FOR  ALL  R  LET  DFCYN3DF2CR) , R) =DELTA3*YN3DF3(R)/A; 

X  SFECIAL  CASE; 

N:  =0; 

X  ASSUMED  SOLUTION  WITH  THE  BESSEL  FUNCTION  OF  SECOND  KIND; 

Wl:=  (AlxJNl(R)  +  A4XYN1(R))XC0SCNXTHETA); 

W2:=  CA2XJN2CR)  +  A5XYN2CR) )*COSCNXTHETA) ; 

W3:=  CA3XJN3CR)  +  A6XYN3C R) JXSINCNXTHETA ) ; 

X  POTENTIAL  FUNCTIONS; 

KSIR==  CSIGMA1-I )XDFCWI,R)  +  CSIGMA2-1)»DF(W2,R)  +  DFC W3, THETA )/R; 
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KSIT:=  (SIGMA1-1 )XDF( HI, THETA )/R  +  <SIGMA2-l)xDF(W2,THETA)/R  - 
DF(H5,R); 

H : -  HI  +  H2 ' 

X*  MOMENTS* AND  SHEARS; 

MR : =  Dx( DF( KSIR, R)  +  (MU/R)X(KSIR  +  DF( KSIT , THETA ))) ; 

MT:=  D*( CKSIR  +  DF( KSIT, THETA) )/R  +  MUXDFCKSIR, R) ) ; 

MRT : =  CDX(1-MU)/2)X( (DF( KSIR, THETA)  -  KSIT)/R  +  DF(KSIT , R) ) ; 

QR : =  K**2*G*H*( KSIR  +  DF(H,R) ) ; 

QT;=  KX*2XGXHX(KSIT  +  DF(M,THETA)/R) ; 

X  BOUNDRY  CONDITIONS; 

X  R  =  A,  MR=0,MRT=0,  AND  QR=0; 

MRBCA : =  SUB(R=A,MR); 

MRTBCA : *  SUB( R=A, MRT) ; 

QRBCA : =  SUB( R=A, QR) ; 

%  R  =  B  AT  THE  CENTER, KSIR=0, KSIT=0,  W  =  HTO; 

X  HTO  IS  THE  DRIVING  POINT  DISPLACEMENT; 

B  s  =  0  ; 

KSIRBCB : =  SUB(R=B,KSIR) ; 

KSITBCB •.  -  SUBC R-B, KSIT ) ; 

WBCB:=  SUB(R=B, W)  -  HTO; 

K  SOLVE  FOR  THE  CONSTANTS  A1  TO  A6; 

AA1 : -  <  -  B4XC34XC23XC12  +  B4xC34»C22xC13  +  B4XC33XC24XC12 

-  B4XC33XC22XC14  -  B4*C32xC24xCl 3  +  B4XC32XC23XC14 
+  C44XB3XC23XC12  -  C44XB3XC22XC13  -  C44XC33XB2XC12 
+  C44XC33XC22XB1  +  C44XC32XB2XC13  -  C44XC32XC23XB1 

-  C43XB3XC24XC12  +  C43XB3XC22XC14  +  C43XC34XB2XC12 

-  C43XC34XC22XB1  -  C43XC32XB2XC14  +  C43XC32XC24XB1 
+  C42XB3XC24XC13  -  C42XB3XC23XCI4  -  C42XC34XB2XC13 

+  C42XC34XC23XB1  +  C42XC33XB2XC14  -  C42XC33XC24XB1 )/ 
(C44XC33XC22XC11  -  C44XC33XC21XC1 2  -  C44XC32XC23XC11  +  C44XC32*C21» 
C13  +  C44XC31XC23XC12  -  C44XC31XC22XC13  -  C43XC34XC22XC1 1  +  C43X 
C34XC21XC12  +  C43XC32XC24XC11  -  C43xC32*C21xC14  -  C43XC31XC24XC12 
+  C43XC31XC22XC14  +  C42XC34XC23XC11  -  C42XC34XC21XC13  -  C42*C33X 
C24XC11  +  C42XC33XC21XC14  +  C42XC31*C24XC13  -  C42XC31XC23XC14  - 
C41XC34XC23XC12  +  C41»C34xC22xC13  +  C41XC33XC24XC12  -  C41XC33XC22X 
C14  -  C41XC32XC24XC13  +  C41XC32XC23XCI4) ; 

AA2:=  (B4XC34XC23XC11  -  B4xC34xC2lxC13  -  B4XC33XC24XC11  +  B4 
XC33XC21XC14  +  B4XC31XC24XC13  -  B4*C31XC23*C14  -  C44X 
B3XC23XC11  +  C44XB3XC21XC13  +  C44XC33XB2XC11  -  C44x 
C33XC21XB1  -  C44XC31XB2XC13  +  C44XC31XC23XB1  +  C43XB3 
XC24XC11  -  C43XB3XC21XC14  -  C43XC34XB2*C1 1  +  C43XC34X 
C21XB1  +  C43XC31XB2XC14  -  C43XC31XC24XB1  -  C41XB3XC24 
XC13  +  C41XB3XC23XC14  +  C41XC34XB2XC13  -  C41XC34XC23X 
B1  -  C41XC33XB2XC14  +  C41XC33XC24XB1 )/CC44xC33xC22x 
Cll  -  C44XC33XC21XC12  -  C44XC32XC23XC11  +  C44XC32XC2IXC13  +  C44X 
C31XC23XC12  -  C44XC31XC22XC13  -  C43XC34XC22XC1 1  +  C43XC34XC21X 
C12  +  C43*C32XC24XC11  -  C43XC32XC21XC14  -  C43XC31XC24XC12  +  C43X 
C31XC22*C14  +  C42XC34XC23XC11  -  C42*C34XC2lxC13  -  C42XC33XC24X 
Cll  +  C42XC33XC21XC14  +  C42XC31XC24XC13  -  C42XC3IXC23XC14  -  C4lx 
C34XC23XC12  +  C41XC34XC22XC13  +  C41XC33XC24XC12  -  C4lxC33*C22X 
C14  -  C41XC32XC24XC13  +  C41XC32XC23XC14) ; 

AA4:=  (  -  B4XC34XC22XC11  +  B4XC34XC21XC12  ♦  B4*C32*C24«C11 

-  B4XC32XC21XC14  -  B4XC31XC24XC12  +  B4XC3IXC22XC14 
+  C44XB3XC22XC11  -  C44XB3XC21XC12  -  C44XC32XB2XC11 
+  C44XC32XC21XB1  +  C44XC31XB2XC12  -  C44*C3lxC22xBl 

-  C42XB3XC24XC11  +  C42XB3XC21XC14  +  C42XC34XB2XC11 

-  C42XC34XC21XBI  -  C42XC31XB2XC14  +  C42XC3IXC24XB1 
+  C4IXB3XC24XC12  -  C41xB3»C22xCI4  -  C41XC34XB2XC12 

+  C41XC34XC22XB1  +  C4IXC32XB2XC14  -  C41XC32XC24XB1 )/ 
(C44XC33XC22XC11  -  C44*C33*C21XC12  -  C44*C32*C23*C11  +  C44*C52*C21* 
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C13  +  C44XC31XC23XC12  -  C44XC31XC22XC13  -  C43XC34XC22XC11  +  C43x 
C34XC21XC12  +  C43XC32XC24XC11  -  C45XC32XC21XC14  -  C43XC31XC24XC12 
+  C43XC31XC22XC14  +  C42XC34XC23XC11  -  C42XC34XC21XC13  -  C42XC33X 
C24XC11  +  C42XC33XC21XC14  +  C42XC31XC24XC13  -  C42XC31XC23XC14  - 
C41XC34XC23XC12  +  C41XC34XC22XC13  +  C41XC33XC24XC12  -  C41XC33XC22X 
C14  -  C41XC32XC24XC13  +  C41XC32XC23XC14 ) ; 

AA5:=  CB4XC33XC22XC11  -  B4XC33XC21XC12  -  B4XC32XC23XC11  +  B4 
XC32XC21XC13  +  B4XC31XC23XC12  -  B4XC31XC22XC13  -  C43X 
B3XC22XC11  +  C43XB3XC21XC12  +  C43XC32XB2XC11  -  C43x 
C32XC21XB1  -  C43XC31XB2XC12  +  C43XC31XC22XB1  +  C42XB3 
XC23XC11  -  C42XB3XC21XC13  -  C42XC33XB2XC11  +  C42XC33X 
C21XB1  +  C42XC31XB2XC13  -  C42XC31XC23XB1  -  C41XB3XC23 
XC12  +  C41XB3XC22XC13  +  C41XC33XB2XC12  -  C41XC33XC22X 
B1  -  C41XC32XB2XC13  +  C41XC32XC23XB1 )/(C44XC33xC22X 
Cll  -  C44XC33XC21XC12  -  C44XC32XC23XC11  +  C44XC32XC21XC13  +  C44X 
C31XC23XC12  -  C44XC31XC22XC13  -  C43XC34XC22XC11  +  C43XC34XC21X 
C 12  +  C43XC32XC24XC11  -  C43XC32XC21XC14  -  C43XC31XC24XC12  +  C43x 
C31XC22XC14  +  C42XC34XC23XC11  -  C42XC34XC21XC13  -  C42XC33XC24X 
Cll  +  C42XC33XC21XC14  +  C42XC31XC24XC13  -  C42XC31XC23XC14  -  C41X 
C34XC23XC12  +  C41XC34XC22XC13  +  C41XC33XC24XC12  -  C41XC33XC22X 
C14  -  C41XC32XC24XC13  +  C41XC32XC23XC14) ; 

A3 i *  AA3; 

A6:=  AA6 ; 

%  DERIVATION  OF  Cll  TO  C44  AND  B1  TO  B4; 

Cll:*  SUB(A1=1,A2=0,A4=0,A5=0,MRBCA); 


C12 : 3  SUB(A1=0,A2=1,A4=0,A5=0,MRBCA) 

Cl 3 : *  SUB(A1=0,A2=0,A4=1,A5=0,MRBCA) 

C14:*  SUB(A1=0,A230,A4=0,A5=1.MRBCA) 

C21:=  SUB(A1=1,A2=0,A4=0,A5=0,0RBCA) 

C22 : 3  SUB(A1=0,A2=1,A4=0,A5=0,QRBCA) 

C23:=  SUB(A1=0,A2=0,A4=1/A5=0, QRBCA ) 

C24:=  SUB (A1=0,A2=0,A4=0,A5S1, QRBCA) i 
C31  •• 3  SUB(A1=1,A230,A4=0,A5=0,KSIRBCB) 
C32:=  SUB(A1=0,A2=1»A4=0,A5=0, KSIRBCB ) 
C33:=  SUB(A1=0,A2=0,A4=1,A5=0, KSIRBCB) 

C34 . 3  SUB(A1=0,A2=0,A4=0,A5=1. KSIRBCB); 

C41 : 3  SUB(A1=1,A2=0,A4=0,A5=0,WTO=0,WBCB>; 
C42:=  SUB(A1=0,A2=1,A4=0,A5=0,WT0=0,WBCB); 
C43:=  SUB(A1=0,A2=0,A4=1,A5=0,WT0=0,HBCB>; 
C44:*  SUB(A1=0,A2=0,A430,A5=1,WT0=0,WBCB); 


Bl:=  0; 

B2 : -  0; 

B3:=  0; 

B4:=  WTO; 

A 1 : 3  AA1 ; 

A2:=  AA2; 

A4 : =  AA4; 

A5 : -  AA5 ; 

F0FT2 : =  INT(RXQR, THETA); 

FOFT1 : =  SUB(R=B,F0FT2); 

FOFT j  =  SUB(THETA  =  2XPI,F0FT1)  -  SUB( THETA3 0 , FOFT1 ) ; 
Z.-=  FOFT/ ( OMEGAXWT 0 )  ; 

Si  SPECIAL  CASE  OF  N  =  0; 

ON  FORT; 

WRITE  Z; 

OFF  FORT; 

SHUT  TT1F$ 


BYE; 
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APPENDIX  B. 

FORTRAN  PROGRAM  TO  CALCULATE  THE  DRIVING  POINT 

IMPEDANCE 


PROGRAM  THESIS 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

cx  THE  PURPOSE  OF  THIS  PROGRAM  IS  TO  GENERATE  THE  DATA  WHICH  IS  x 

CX  NECESSARY  TO  CALCULATE  THE  DRIVING  POINT  IMPEDANCE  VS  FREQUENCY  x 

CXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


Cx 

LIST  OF 

VARIABLES 

X 

Cx 

FL 

LOWER  FREQUENCY  OF  INTEREST,  HZ 

X 

CX 

FH 

UPPER  FREQUNECY  OF  INTEREST,  HZ 

X 

CX 

F 

FREQUENCY  OF  INTEREST,  HZ 

X 

Cx 

OMEGA 

ANGULAR  FREQUENCY  OF  INTEREST,  =2PIXF,  RAD/SEC 

X 

CX 

A 

RADIUS  OF  OUTER  EDGE  OF  PLATE,  IN 

X 

Cx 

B 

RADIUS  OF  WAVE  GUIDE  ATTACHMENT,  IN 

X 

CX 

H 

THICKNESS  OF  PLATE,  IN 

X 

Cx 

TEMP 

AVERAGE  PLATE  TEMP  IN  DEGREES  F 

X 

CX 

IMAT 

TYPE  OF  MATERIAL 

X 

Cx 

=  0,  ELASTIC 

X 

CX 

=  1,  VISCOELASTIC  WITH  DATA  FROM  G.G.LEE 

X 

Cx 

=  2,  VISCOELASTIC  WITH  DATA  FROM  UNITED  MCGILL 

X 

CX 

RHO 

DENSITY  OF  PLATE,  SLUGS/INXX3 

X 

Cx 

E 

YOUNG'S  MODULUS,  PSI 

X 

CX 

G 

SHEAR  MODULUS, PSI 

X 

CX 

MU 

POISSON'S  RATIO 

X 

cx 

JXY 

BESSEL  FUNCTION  OF  THE  FIRST  KIND,  EVALUATED  AT 

X 

CX 

DELTAXXY/A 

X 

Cx 

YXY 

BESSEL  FUNCTION  OF  THE  SECOND  KIND,  EVALUATED  AT 

X 

CX 

DELTAXXY/A 

X 

Cx 

JXDF1Y 

FIRST  DERIVATIVE  OF  THE  BESSEL  FUNCTION  OF  THE 

X 

Cx 

FIRST  KIND,  EVALUATED  AT  DELTAXXY/A 

X 

Cx 

YXDF1Y 

FIRST  DERIVATIVE  OF  THE  BESSEL  FUNCTION  OF  THE 

X 

CX 

SECOND  KIND,  EVALUATED  AT  DELTAXXY/A 

X 

CX 

JXDF2Y 

SECOND  DERIVATIVE  OF  THE  BESSEL  FUNCTION  OF  THE 

X 

Cx 

FIRST  KIND,  EVALUATED  AT  DELTAXXY/A 

X 

CX 

YXDF2Y 

SECOND  DERIVATIVE  OF  THE  BESSEL  FUNCTION  OF  THE 

X 

CX 

SECOND  KIND,  EVALUATED  AT  DELTAXXY/A 

X 

CX 

DELTAI 

DEMENSIONLESS  PARAMETER 

X 

Cx 

DELTA2 

DEMENSIONLESS  PARAMETER 

X 

CX 

DELTA3 

DEMENSIONLESS  PARAMETER 

X 

Cx 

SIGMA1 

DEMENSIONLESS  PARAMETER 

X 

CX 

SIGMA 2 

DEMENSIONLESS  PARAMETER 

X 

Cx 

DRIVING  POINT  IMPEDANCE 

X 

Cx 

ZZ 

COMPLEX  IMPEDANCE,  LBF-SEC/IN 

X 

CX 

z 

MAGNITUDE  OF  ZZ,  LBF-SEC/IN 

X 

Cx 

ZR 

REAL  PART  OF  ZZ,  LBF-SEC/IN 

X 

CX 

ZI 

IMAGINARY  PART  OF  ZZ,  LBF-SEC/IN 

X 

CX 

ZDB 

Z  IN  DB, 

X 

CX 

K 

SHEAR  COEFFICEINT 

X 

Cx 

PI 

3.14 - 

X 

C* 

X 

DUMMY  ARGUMENT  FOR  BESSEL  FUNCTION  EVALUATION 

X 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXX 

IMPLICIT  C0MPLEXX8  (A,D,J,S,Y) 
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on  ooo  o no  nnnn 


COMPLEXXS  E,X,G 
REAL  MU,A,K 

CHARACTER  ANS*1,  NAMEX6 

C0MM0N/BESS1A/J1A, J1DF1A, J1DF2A, Y1A. Y1DF1A»Y1DF2A 
C0MM0N/BESS1B/J1B, J1DF1B, J1DF2B, Y1B, Y1DF1B.Y1DF2B 
C0MM0N/BESS2A/J2A, J2DF1A, J2DF2A, Y2A, Y2DF1A, Y2DF2A 
C0MM0N/BESS2B/J2B, J2DF1B, J2DF2B , Y2B, Y2DF1B, Y2DF2B 
C0MMQN/BESS3A/J3A, J3DF1A, J3DF2A, Y3A, Y3DF1A, Y3DF2A 
C0MM0N/BESS3B/ J3B, J3DF1B, J3DF2B, Y3B, Y3DF1B, Y3DF2B 
COMMON/GEOMTY/ A , B , H , K 
COMMON/MAT/ IMAT , E , G , MU , RHO ,  T EMP 
COMMON/ CONST 1/ DEL TA1 ,  DELTA2, DELTA3, SIGMA1 , SIGMA2 
PI=  ARCOS(-l.O) 

WRITE!6,X)'  PI  =  ',PI 

INPUT  OF  MATERIAL  TYPE 


5 


CONTINUE 

CALL  FRTCMS  ! 'CLRSCRN  ') 

WRITE!6,x)»  THE  FOLLOWING  MATERIALS  ARE  SUPPORTED' 

WRITEC  6  X ) f  1 

WRITE(6/*) '  TYPE  NAME  CHARACTERISTICS  FROM 

1* 


WRITEC 6  X ) *  ' 

WRITEC6»X) 1  ELASTIC! E)  ALUMINUM  GERE-TIMOSHENKO • 

WRITEC 6  x ) •  ' 

WRITEC6^X) •  VISCOELASTIC! V)  LD-400  VIBRATION  DAMPING' 

WRITEC6 , X ) '  VISCOELASTIC! M)  C-2204  UNITED  MCGILL  • 

WRITE! 6  X ) *  ' 

WRITEC6ix)'  WHAT  TYPE  OF  MATERIAL  IS  UNDER  STUDY?' 

WRITE! 6 , ») '  TYPE  E,  V,  OR  M.* 

READ! 5. ' ( A1 ) • )  ANS 
WRITE! 6, X)'  ' , ANS 

IF  !  ANS  .HE.  'E*  .AND.  ANS  .NE.  'V'  .AND.  ANS  .NE.  'M*  )  GO  TO  5 
IMAT  =  0 

IF  !  ANS  .EQ.  »V»  )  IMAT  *  1 
IF  C  ANS  .EQ.  'M'  )  IMAT  =  2 
WRITEC 6 , X) '  IMAT= ' , IMAT 


INPUT  OF  PLATE  TEPERATURE 

6  WRITE! 6 , X) »  WHAT  IS  THE  AVERAGE  TEMPERATURE  OF  THE  PLATE  IN  DEGREE 
1  F7CTEMP) ' 

READ! 5 ,  * )  TEMP 

WRITE!6,X) '  TEMP= ' , TEMP, '  F' 

IF  C  TEMP  .LT.  63.0  .OR.  TEMP  .GT.  104.0  )  GO  TO  6 


INPUT  OF  PLATE  GEOMETRY 

WRITE!6,X) '  WHAT  IS  THE  RADIUS  OF  THE  PLATE  IN  INCHSTCA)' 
READC5.X)  A 

WRITEC6.X) '  A= ' , A, '  IN' 

IF  C  A  .LE.  0.0  )  GO  TO  5 
CONTINUE 

WRITEC 6 , X ) '  WHAT  IS  THE  MOUNT  RADIUS  IN  INCHSTCB)' 

READ! 5,X)  B 
B  =  0.3125 

WRITE! 6  *  * ) '  B=’,B» '  IN' 

IF  !  B  .LE.  0.0  )  GO  TO  7 
8  CONTINUE 
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WRITE(6,X) •  WHAT  IS  THE  THICKNESS  OF  THE  PLATE  IN  INCHS?(H)' 
READ(5,x)  H 
H  =  3. 0/8.0 

WRITE(6,x) *  H=',H, •  IN* 

IF  (  H  .LE.  0.0  )  GO  TO  8 

INPUT  OF  FREQUENCY  RANGE 

10  WRITE(6,X)'  WHAT  IS  THE  LOW  FREQUENCY  OF  INTEREST,  IN  HZ?(FL)' 
READ(5,x)  FL 

WRITE(6 , x) 1  FL=*,FL, '  HZ' 

IF  (  FL  .LE.  0.0  )  GO  TO  10 

15  WRITE(6,x) '  WHAT  IS  THE  HIGH  FREQUENCY  OF  INTEREST,  IN  HZ?(FH) ' 
READ(5, *)  FH 

WRIT£(6,x)'  FH= ' , FH, •  HZ' 

IF  (  FH  .LE.  FL  >  GO  TO  15 

INPUT  OF  LOCATION  TO  WRITE  RESULTS 

WRITE(6 , X)  '  WHAT  IS  THE  FILE  NAME?  ' 

READC6, ’(A)')  NAME 

OPEN ( UNIT3 15, FILE3  NAME,  F0RM=  'FORMATTED' , ACCESS3  'SEQUENTIAL ', ST 
1ATUS=  'NEW'  ) 

IF  (  IMAT  .EQ.  0  )  IP0INT  =  99 
IF  C  IMAT  .EQ.  1  .OR.  IMAT  .EQ.  2  )  IPOINT  =  99 

START  OF  CALCULATIONS 

DO  90  1=0, IPOINT 

FREQUENCY  OF  INTEREST 

F  3  FL  +  (FH  -  FL ) XFLOAT ( I )/ FLOAT (IPOINT ) 

OMEGA  =  2 . 0XPI#F 

WRITEC6 , X) '  I=',I,»  F=',F,'HZ  OMEGA3 ', OMEGA, 'RAD/SEC* 

GENERATION  OF  DIMENSIONLESS  PARAMETERS 

CALL  CONST(OMEGA, F) 

GENERATION  OF  BESSEL  FUNCTIONS 

BESSEL  FUNCTIONS  FOR  DELTA1  AND  RADIUS  A 
X  =  DELTA1 
C  WRITE(6 , X ) '  x= ' , X 

CALL  BESSEL (X. J1 A, JlDFlA, J1DF2A, Y1A , Y1 DF1A, Y1DF2A) 

C  BESSEL  FUNCTIONS  FOR  DELTA1  AND  RADIUS  B 
X  =  DELTA1*B/A 
C  WRITE(6,X) '  X=',X 

CALL  BESSEL (X, JIB , J1DF1B , J1DF2B, Y1B , Y1DF1B, Y1DF2B) 

C  BESSEL  FUNCTIONS  FOR  DELTA2  AND  RADIUS  A 
X  =  DELTA2 

C  WRITE(6 , X) •  X=*,X 

CALL  BESSEL(X, J2A, J2DF1A, J2DF2A, Y2A, Y2DF1A,Y2DF2A) 

C  BESSEL  FUNCTIONS  FOR  DELTA2  AND  RADIUS  B 
X  =  DELTA2XB/A 
C  WRITE(6 , X) »  X=',X 

CALL  BESSEL(X, J2B, J2DF1B, J2DF2B, Y2B, Y2DF1B,Y2DF2B) 

C  BESSEL  FUNCTIONS  FOR  DELTA3  AND  RADIUS  A 

X  =  DELTA3 
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C  WRITEC6 * *) 1  Xs*  X 

CALL  BESSEL<X,J3A,J3DF1A,J3DF2A,Y3A,Y3DF1A,Y3DF2A) 

C  BESSEL  FUNCTIONS  FOR  DELTA3  AND  RADIUS  B 

X  =  DELTA3KB/A 
C  WRITEC6**) 1  x= 1 9 X 

CALL  BESSEL (X,  J3B, J3DF1B, J3DF2B,Y3B,Y3DF1B,Y3DF2B) 

C 

C  IMPERANCE  CALCULATIONS 

C 

CALL  IMPEDOCOMEGA,Z,ZR,ZI,ZDB) 

WRITEC 15,20)  F,Z,ZR,ZI,ZDB 
90  CONTINUE 

WRITE(6,x) '  DO  YOU  WANT  TO  CONTINUETCY/N) . » 

READL5, '(A)')  ANS 
CALL  FRTCMS  ( 'CLRSCRN  •) 

IF  (  ANS  .EQ.  'Y*  )  GO  TO  5 
20  F0RMAT<1X,5(G15.7,1X)) 

END 

SUBROUTINE  CONSTCOMEGA, F) 

CXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXKXXXXXXXXXXXXXXX 


cx 

THE  PUPOSE 

OF  THIS  PROGRAM  IS  TO  DETERMINE  THE  DIMENSIONLESS 

X 

Cx 

PARAMETERS 

AND  YOUNG’S  AND  SHEAR  MODULUS  FOR  THE 

X 

Cx 

VISCOELASTIC  MATERIAL 

X 

CxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxXxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

CX 

LIST  OF  VARIABLE 

X 

cx 

A 

PLATE  RADIUS,  IN 

X 

cx 

B 

WAVEGUIDE  ATTACHMENT  RADIUS,  IN 

X 

cx 

H 

PLATE  THICKNESS,  IN 

X 

cx 

TEMP 

AVERAGE  PLATE  TEMP  IN  DEGREES  F 

X 

cx 

IMAT 

TYPE  OF  MATERIAL 

X 

cx 

=  0,  ELASTIC  MATERIAL 

X 

cx 

=  1.  VISCOELASTIC  WITH  DATA  FROM  G.G.LEE 

X 

cx 

=  2,  VISCOELASTIC  WITH  DATA  FROM  UNITED  MCGILL 

X 

cx 

E 

YOUNG'S  MODULUS,  PSI 

X 

CX 

G 

SHEAR  MODULUS,  PSI 

X 

CX 

ETHA 

VISCOELASTIC  LOSS  FACTOR 

X 

cx 

GMAG 

VISCOELASTIC  SHEAR  MODULUS 

X 

cx 

D 

FLEXURAL  RIGIDITY,  PSI-IN 

X 

CX 

MU 

POISSON'S  RATIO 

X 

CX 

RHO 

MASS  DENSITY  OF  PLATE,  SLUGS/INXX3 

X 

CX 

K 

SHEAR  COEFFICIENT 

X 

CX 

OMEGA 

CIRCULAR  FREQUENCY  OF  INTEREST,  RAD/SEC 

X 

cx 

F 

FREQUENCY  OF  INTEREST,  HZ 

X 

cx 

DELTA1 

DEMENSIONLESS  PARAMETER 

X 

CX 

DELTA2 

DEMENSIONLESS  PARAMETER 

X 

CX 

DELTA3 

DEMENSIONLESS  PARAMETER 

X 

cx 

SIGMA1 

DEMENSIONLESS  PARAMETER 

X 

cx 

SIGMA 2 

DEMENSIONLESS  PARAMETER 

X 

cx 

DELT12 

DELTA1XX2 

X 

CX 

DELT22 

DELTA2XX2 

X 

CX 

DELT32 

DELTA3XX2 

X 

cx 

RR 

DEMENSIONLESS  PARAMETER 

X 

cx 

S 

DEMENSIONLESS  PARAMETER 

X 

cx 

L4 

DEMENSIONLESS  PARAMETER 

X 

cx 

PI 

3.14 _ 

X 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

C0MPLEXX8  DELTA1 , DELTA2, DELTA3, SIGMA1 , SIGMA2, L4, D, E,G, DELT12, 
1DELT22, DELT32 
REAL  MU,  K 

COMMON/GEOMTY/A , B , H , K 


COMMON/MAT/ IMAT , E , G, MU , RHO , TEMP 

COMMON/ CONST1/ DEL TAX , DELTA2, DELTA3, SIGMA1 , SIGMA2 

PIsACOS!-l .0) 

C  WRITE(X,X)'  PI  =  ',PI 

C 

C  CALCULATION  OF  THE  MATERIAL  CHARACTERISTICS 

C 

CALL  YOUNGS! F) 

C 

C  CALCULATION  OF  FLEXURAL  RIGIDITY 

C 

D  =  E*HXX3/<12.0*(1.0-MUX*2)) 

C 

C  SHEAR  COEFFICIENT 

C 

C  K  =  SQRT(PIXX2/12.0) 

C  K  =  0.95541 

K  =  0.92736 

C  WRITE!6 , x) *  K=',K 

C  WRITE!6 , X) *  D* ' ,  D 

C  WRITEC6,x) •  G=',G 

C  WRITE!6,x) '  RHOSRHO 

C 

C  DIMENSIONLESS  PARAMETERS 

C 

RR  =  CH/A)XX2/12.0 
C  WRITE(6 , X ) *  RR= • , RR 

S  =  (1.0/(6.0XKXX2X(1.0-MU)))X(H/A)X*2 
C  WRITE! 6 , X) »  S=',S 

L4  =  RHOXHXAXK4XQMEGAXX2/!12.0XD) 

C  WRITE! 6 , X) '  L4=',L4 

DELT12  =  !L4X(RR+S+SQRT!!RR-S)XX2+4.0/L4))/2.0) 

C  WRITE! 6 , x) •  DELT12= ' , DELT12 

DELT22  =  !L4x(RR+S-SQRT! 'RR-S)XX2+4.0/L4))/2.0) 

C  WRITE!6,X) •  DELT22= ' , DELT22 

DELT32  =  !2.0X(RR*L4-1 . 0/S)/ ! 1 . O-MU) ) 

C  WRITE! 6 , X) •  DELT32= ' , DELT32 

SIGMA1  =  DELT22  /! RRXL4-1 . 0/S) 

C  WRITE!6, x) •  SIGMA1 = ' , SIGMA1 

SIGMA2  =  DELT12/!RRXL4-1.0/S) 

C  WRITE! 6 , x ) •  S I GMA  2  =  ' ,  S I GMA  2 

DELTA1  =  SQRT ( DELT12 ) 

C  WRITE! 6 , x ) •  DELTA1= ' , DELTA1 

DELTA2  =  SQRTCDELT22) 

C  WRITE! 6 , X ) ’  DELTA25 ' ,  DELTA2 

DELTA3  =  S0RTCDEIT32) 

C  WR I T  E ! 6 , * ) '  DELTA3=  * , DELTA3 

C  GO  TO  10 

999  CONTINUE 
RETURN 
END 

SUBROUTINE  YOUNGS! F) 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXKKXXXXKKXKKKXKKKKXKXKXXXX 

cx  THE  MATERIAL  CHARACTERISTICS 

CX  THE  PURPOSE  OF  THIS  PROGRAM  IS  TO  GENERATE  YOUNG'S  MODULUS 
CX  AND  SHEAR  MODULUS  FOR  THE  MATERIAL  UNDER  STUDY 
CX  C.D.HETTEMA  4  AUG  88 

CXXXXXXXXXXXXXXXXXKXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX! 

cx  LIST  OF  VARIABLES 

CX  F  FREQUENCY  OF  INTEREST,  HZ 
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xixvximii 


Cx 

FLOG 

LOGIO(F) 

X 

Cx 

TEMP 

AVERAGE  PLATE  TEMP  IN  DEGREES  F 

X 

CX 

IMAT 

TYPE  OF  MATERIAL 

X 

Cx 

=  0,  ELASTIC 

X 

Cx 

=  1,  VISCOELASTIC  WITH  DATA  FROM  G.G.LEE 

X 

Cx 

*  2,  VISCOELASTIC  WITH  DATA  FROM  UNITED 

MCGILL  x 

CX 

RHO 

DENSITY  OF  PLATE,  SLUGS/INXX3 

X 

Cx 

E 

YOUNG’S  MODULUS,  PSI 

X 

Cx 

G 

SHEAR  MODULUS, PSI 

X 

Cx 

ETHA 

LOSS  FACTOR 

X 

Cx 

EP 

REAL  PART  OF  COMPLEX  YOUNG’S  MODULUS 

X 

Cx 

EDP 

COMPLEX  PART  OF  COMPLEX  YOUNG’S  MODULUS 

X 

Cx 

EPX 

REAL  PART  OF  COMPLEX  YOUNG’S  MODULUS  AT 

TEMP  X  x 

Cx 

EDPX 

COMPLEX  PART  OF  COMPLEX  YOUNG’S  MODULUS 

AT  TEMP  X  X 

Cx 

DEP 

DUMMY  VARIABLE  FOR  EP 

X 

CX 

DEDP 

DUMMY  VARIABLE  FOR  EDP 

X 

Cx 

DEPX 

DUMMY  VARIABLE  FOR  EPX 

X 

Cx 

DEDPX 

DUMMY  VARIABLE  FOR  EDPX 

X 

Cx 

MU 

POISSON'S  RATIO 

X 

Cx 

PI 

3.14 - 

X 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXX 
COMPLEXXS  E,G 
REAL  MU 

CHARACTER  ANSX1 

COMMON/MAT/ IMAT , E , G , MU , RHO , TEMP 
C 

C  FOR  USE  IN  UNITED  MCGILL  DATA 

C 

FLOG  =  LOGIO(F) 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXK 

c 

C  ELASTIC  MATERIAL  DATA 

c 

IF  (  IMAT  .EQ.  0  )  THEN 
E  =  (11.0E6,0.0) 

MU  *  0.33 

G  =  E/(2.0X(1.0+MU)) 

RHO  =  3. 06713E-3 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXKXXXXXXXXXXKXXXX 

c 

C  VISCOELASTIC  MATERIAL  CHARACTERISTICS  FROM  NASHIF,  JONES,  AND 

C  HENDERSON 

C 

ELSE  IF  C  IMAT  .EQ.  1  )  THEN 

ETHA  =  0.65xEXP(-0.52732x((ABSCLOG10CF/62)))xxl.956)) 

GMAG  =  0.00002503XFXX3  -  0.1752XFXX2  +  457.5883*F  +  29280 
G  =  GMAGX( ( 1 . 0, 0 . 0 )  +  (0.1,1. 0  JXETHA) 

MU  =  0.5 

E  =  GX(2 . 0*( 1 . 0+MU) ) 

RHO  =  1.71101E-3 

CXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXX 

c 

C  VISCOELASTIC  MATERIAL  CHARACTERISTICS  FROM  UNITED  MCGILL 
C 

ELSE  IF  (  IMAT  .EQ.  2  )  THEN 
C 

C  FOUR  POINT,  THIRD  ORDER  CURVE  FITS  FOR  A  SET  TEMP 

C 

C 

C  TEMP  OF  68  F 
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c 

c 


c 


c 

c 

c 

c 


c 


c 

c 

c 

c 


c 


c 

c 

c 


c 

c 

c 


DEP  FOR  T=68F 
A3*-0 . 129999657E-01 
A2=  0.999998950E-01 
A1=-0.625001788E-01 
AO=  2 . 02499962 

DEDPPFOR  T-68FA1*FL0G  +  A2*FL0G**2  +  A3XFLOGXX3 
A3=-0. 166667 290 E-01 
A2=  0.112500608 
A1=-0.708351135E-01 
A0=  1.77500057 

DEDPA  =  AO  +  AlxFLOG  +  A2XFL0GXX2  +  A3XFL0GXX3 
TEMP  OF  86  F 


DEP  FOR  T=86F 
A3=  0.250002518E-01 
A2=-0. 187501898 
Al=  0.612509189 
A0=  0.999996593 

DEDPPFOR  T-86FA1XFL0G  +  A2*FL0G**2  +  A3*FL0G*X3 
A3=  0.129999657E-01 
A2=-0.999993682E-01 
Al=  0.387997609 
A0=  1.05000119 

DEDPB  =  AO  +  AlxFLOG  +  A2XFL0GXX2  +  A3XFL0GXX3 
TEMP  OF  109  F 


DEP  FOR  T*1 09F 

A3=  0.158995682E-06 
A2=  0 . 129988556 E-01 
Al=  0.87 5023007 E-01 
A0=  0.899998665 

DEDPPFOR  f=109F1XFL°G  +  A2*FL0G**2  +  A3XFL0GXX3 
A3=  0.833320618E-02 
A2=-0. 62999096 3 E-01 
Al=  0.259169696 
A0=  0.800001199 

DEDPC  =  AO  +  AlxFLOG  +  A2XFL0GXX2  +  A3xFL0Gxx3 
LINEAR  FIT  FOR  SET  FREQUENCY 


THEN 


IF  (  TEMP  .GE.  68.0  .AND.  TEMP  .LT.  86.0  ) 

DEP  =  ( DEPB-DEPA  )  X(  TEMP-68 .0/18.0  +  DEPA 
DEDP  =  (DEDPB- DEDPA) x(TEMP-6 8 .01/18.0  +  DEDPA 
n«E  lF,  LTEMP  ,GE-  86 -°  -AND-  TEMP  -LE.  109.0 
DEP  *  (DEPC-DEPB)X(T£MP-86.0)/18.0  +  DEPB 

EHDPIF  <DEDPC~DEDPB)*(TEMP_86 . 0V1S .0  +  DEDPB 


CONVERSION  OF  DATA  TO  PSI  UNITS 

EP  *  1 .0E8X10XXDEPX1 .95097E-5 
EDP  =  1.QE8X10XXDEDPX1.95097E-5 

ETHA^s^EDP/EP  *  “P«<0.0.1.0, 


THEN 
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MU  =  0.5 

G  =  E/( 2 . 0*( 1 . O+MU) ) 

RHO  »  1 . 70736E-3 

END  IF 

RETURN 

END 

SUBROUTINE  BESSEL (XX, JO, JODF1,JODF2,YO.YODF1,YODF2) 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


Cx 

THIS  PROGRAM 

CALCULATES  THE  BESSEL  FUNCTION  OF  THE  FIRST  AND 

X 

CX 

SECOND 

KIND  , 

ALONG  WITH  THE  FIRST  AND  SECOND  DERIVATIVE  OF  EACH 

X 

Cx 

THIS  WILL  BE 

BY  TAYLOR  SERIES  EXPANTION 

X 

Cx 

N  =  0 

X 

CXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXKXXXXXXXXXX 

CX 

LIST  OF 

VARIABLES 

X 

cx 

M 

COUNTER 

X 

CX 

Ml 

COUNTER  FOR  DERIVATIVES,  M1=M+1 

X 

CX 

ZM 

DOUBLE  PRECISION  M 

X 

CX 

ZM1 

DOUBLE  PRECISION  Ml 

X 

CX 

ZMI 

SUM  OF  THE  INVERSE  OF  M 

X 

cx 

X 

ARGUMENT  OF  THE  FUNCTION 

X 

CX 

XXXD 

DOUBLE  PRECISION  VALUE  OF  XXX 

X 

CX 

JO 

BESSEL  FUNCTION,  FIRST  KIND 

X 

cx 

J0DF1 

FIRST  DERIVATIVE  OF  JO 

X 

cx 

J0DF2 

SECOND  DERIVATIVE  OF  JO 

X 

cx 

JOA 

CHANGE  IN  JO  EACH  STEP 

X 

cx 

JDF1A 

CHANGE  IN  J0DF1  EACH  STEP 

X 

cx 

JDF2A 

CHANGE  IN  J0DF2  EACH  STEP 

X 

CX 

J01 

ABS  OF  JOA 

X 

CX 

JDF11 

ABS  OF  JDF1A 

X 

cx 

JDF2I 

ABS  OF  JDF2A 

X 

CX 

YO 

BESSEL  FUNCTION,  SECOND  KIND 

X 

cx 

Y0DF1 

FIRST  DERIVATIVE  OF  YO 

X 

cx 

Y0DF2 

SECOND  DERIVATIVE  OF  YO 

X 

cx 

YOA 

CHANGE  IN  YO  EACH  STEP 

X 

cx 

YDF1A 

CHANGE  IN  YODFI  EACH  STEP 

X 

CX 

YDF2A 

CHANGE  IN  Y00F2  EACH  STEP 

X 

CX 

Y01 

ABS  OF  YOA 

X 

cx 

YDF11 

ABS  OF  YDF1A 

X 

CX 

YDF21 

ABS  OF  YDF2A 

X 

CX 

GAMMA 

EULER'S  CONSTANT 

X 

cx 

ERRORM 

MAXIMUM  ALLOWED  ERROR 

X 

cx 

ZZZMAX 

MAXIMUM  SIZE  NUMBER  ALLOWED 

X 

cx 

MCI , MC2 

,  MC3 

COUNTERS 

X 

cx 

D 

DUMMY  VARIABLE 

X 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

c 

C  QUAD  COMPLEX 

c 

IMPLICIT  REALX16  (A-H,0-Z) 

c 

C  DOUBLE  COMPLEX 

C 

C  IMPLICIT  C0MPLEXX16  (A-H,0-Z) 

C0MPLEXX32  JOD, JODFID, J0DF2D, JOAD, JDF1AD, JDF2AD, J01D, JDF11D, JDF21d 
C0MPLEXX32  YOD, YODF1D, Y0DF2D, YOAD, YDF1AD, YDF2AD, Y01D, YDF11D, YDF21D 
COMPLEXX32  X 

C0MPLEXX8  JO, JODF1, J0DF2 
COMPLEXxg  Y0,Y0DF1,Y0DF2 
C0MPLEXX8  XX 
CHARACTER  ANSXl 
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non  non  non  non  non  noon  non  non 


ZERO  =  0.0D0 
ONE  =  1.0D0 
TWO  =  2.0DO 
PI  »  DARCOS(-ONE) 

MAXIMUM  ERROR  ALLOWED 

ERRORM  =  1.0D-15 

MAXIMUM  SIZE  NUMBER  ALLOWED 

ZZZMAX  =  1.0D30 

GAMMA  =  0.57721 56 64901 53286 06 06 512D0 

CONTINUE 

JO  ~  ZERO 

CONVERSION  OF  X  TO  QUAD  PERCISION 
X  =  XX  *  ONE 

CHECK  TO  PREVENT  YO  FROM  BEING  UNBOUNDED 
D  =  ABS(X) 

IF  (  D  .LT.  1.0D-3  )  THEN 

WRITEL6 , X) ’  ****  X  IS  TOO  SMALL  ***** 

GO  TO  50 
END  IF 
M  =  0 
CONTINUE 

JO  CALCULATIONS 

Ml  =  M  +  1 
ZM  =  DBLE(M) 

ZM1  =  DBLECM1) 

SET  UP  OF  TAYLOR  SERIES 

IF  (  M  ,EQ.  0  )  THEN 

J0D  =  ONE 

J0AD  =  ZERO 

J01D  =  ONE 

J0DF1D  =  -ONEXX/TWO 

JDF1AD  =  ZERO 

JDF11D  =  ONEXX/TWO 

J0DF2D  =  -ONE/TWO 

JDF2AD  =  ZERO 

JDF21D  =  ONE/TWO 

MLAST  =  0 

MCI  =  1 

MC2  =  1 

MC3  =  1 

EXPANSION  OF  TAYLOR  SERIES 

ELSE  IF  <  M  .GT.  0  )  THEN 
M2  =  M/2 

FOR  JO 


122 


non  oon  non  non  non 


DA  -  ABSCJ01D) 

IF  (  D  .GT.  ERRORM  )  THEN 
J01D  =  J01DX((X/(TWOXZM))XX2) 

JOAD  =  J01D 

IF  (  M2  .EQ.  MLAST  )  JOAD  =  -J01D 
JOD  =  JOD  +  JOAD 
MCI  =  MCI  +  1 
END  IF 

FOR  JODF 

DB  =  ABS( JDF11D) 

IF  (  D  .GT.  ERRORM  )  THEN 

JDF11D  =  JDF11DXXXX2/(CTW0X(ZM1-0NE))X(TW0XZM1)) 

JDF1AD  =  -  JDF11D 

IF  (  M2  .EQ.  MLAST  )  JDF1AD  =  JDF1ID 
J0DF1D  =  J0DF1D  +  JDF1AD 
MC2  -  MC2  +  1 
END  IF 

FOR  J0DF2 

DC  =  ABSC  JDF21D) 

IF  (  D  .GT.  ERRORM  )  THEN 

JDF21D  =  JDF21D*XX*2*(TWO*ZM1-QNE)/((TWO*(ZM1-ONE>)X(TWO*ZM1 
l)XCTWOX(ZMl-ONE)  -  ONE)) 

JDF2AD  =  -  JDF21D 

IF  (  M2  .EQ.  MLAST  )  JDF2AD  =  JDF21D 

J0DF2D  =  J0DF2D  +  JDF2AD 

MC3  s  MC3  +  1 

END  IF 

MLAST  =  M2 

END  IF 

CHECK  IF  SERIES  CONVERGED 

D  =  ABS(JOD) 

D1  =  ABS(JODFID) 

D2  =  ABSC  J0DF2D) 

IF  (  M  .GT.  300  .OR.  D  .GT.  ZZZMAX  .OR.  DA  .GT.  ZZZMAX  .OR. 
1D1  .GT.  ZZZMAX  .OR.  DB  .GT.  ZZZMAX  .OR. 

2D2  .GT.  ZZZMAX  .OR.  DC  .GT.  ZZZMAX  )  THEN 
HRITEC6 , X) '  xxxx  FUNCTION  HAS  NOT  CONVERGED  xxxx* 

GO  TO  50 
END  IF 

IF  C  DA  .LE.  ERRORM  .AND.  M  .GT.  0  .AND. 

1DB  .LE.  ERRORM  .AND.  DC  .LE.  ERRORM  )  GO  TO  20 
M  =  M  +  1 
GO  TO  60 
CONTINUE 

CONVERSION  FROM  QUAD  TO  SINGLE  PERCISION 

JO  =  J0DX1.0 
J0DF1  =  J0DF1DX1.0 
J0DF2  =  J0DF2DX1.0 

YO  CALCULATIONS 

M  *  0 
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oon  ooo  noo  ooo  ooo  ooo 


SET  UP  OF  TAYLOR  SERIES 

YOD  =  (TWO/PI )*<LOG(X/TWO)+GAMMA)*JOD 

Y0DF1D  =  CTWO/PI)X((IOG<X/TWO)+GAMMA)XJODF1D  +  JOD/X)  +  X/PI 
Y0DF2D  =  (TWO/PI )X((L0G(X/TW0)+GAMMA)XJ0DF2D  +  TWOXJODFID/X  - 
1J0/(XXX2))  +  ONE/PI 
Y01D  =  ONE 
YDF11D  =  X/TWO 
YDF21D  =  ONE/TWO 
ZMI  =  ZERO 
MLAST  =  0 
MCI  =  0 
MC2  =  1 
MC3  =  1 

EXPANSION  OF  TAYLOR  SERIES 

30  M  =  M  +  1 

ZM  =  DBLE(M) 

ZMI  =  ZMI  +  ONE/ZM 
M2  =  M/2 

FOR  Y0 

DA  =  ABS(YOID) 

IF  (  DA  .GT.  ERRORM  )  THEN 
Y01D  =  YQlDX((X/(TW0xZMmx2) 

YOAD  =  -Y01D 

IF  (  M2  .EQ.  MLAST  )  YOAD  =  Y01D 
YOD  =  YOD  +  TWOXYOADXZMI/PI 
MCI  =  MCI  +  1 
END  IF 

FOR  YODF 

DB  =  ABS(YDFllD) 

IF  (  DB  .GT.  ERRORM  .AND.  M  .GT.  1  )  THEN 
YDF11D  =  YDF11DXXX*2/((TW0X(ZM-0NE))X(TW0XZM)) 

YDF1AD  =  -YDF11D 

IF  (  M2  .EQ.  MLAST  )  YDF1AD  =  YDF11D 
Y0DF1D  =  Y0DF1D  +  YDF1ADXTN0XZMI/PI 
MC2  =  MC2  +  I 
END  IF 

FOR  Y0DF2 

DC  =  ABS( YDF21D) 

IF  (  DC  .GT.  ERRORM  .AND.  M  .GT.  1  )  THEN 

YDF21D  =  YDF21DXXX*ZXCTWOXZM-ONE)/(CTWOX(ZM-ONE))X(TWOXZM) 
lx(TWOXCZM-ONE)  -  ONE)) 

YDF2AD  =  -YDF2ID 

IF  (  M2  .EQ.  MLAST  )  YDF2AD  =  YDF21D 
Y0DF2D  =  Y0DF2D  +  YDF2ADXTW0XZMI/PI 
MC3  =  MC3  +  1 
END  IF 

CHECK  IF  SERIES  CONVERGED 
MLAST  =  M2 
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D  =  ABS(YOD) 

D1  =  ABS(YODFID) 

D2  =  ABSCYQDF2D) 

IF  (  M  .GT.  300  .OR.  D  .GT.  ZZZMAX  .OR.  DA  .GT.  ZZZMAX  .OR. 

1D1  .GT.  ZZZMAX  .OR.  DB  .GT.  ZZZMAX  .OR. 

2D2  .GT.  ZZZMAX  .OR.  DC  .GT.  ZZZMAX  )  THEN 
WRITE(6,X) *  xxxx  FUNCTION  HAS  NOT  CONVERGED  XX**' 

GO  TO  40 
END  IF 

IF  (  DA  .LE.  ERRORM  .AND.  M  .GT.  0  .AND. 

1DB  .LE.  ERRORM  .AND.  DC  .LE.  ERRORM  )  GO  TO  40 
GO  TO  30 
40  CONTINUE 
C 

C  CONVERSION  FROM  QUAD  TO  SINGLE  PERCISION 
C 

YO  =  Y0DX1.0 
Y0DF1  =  Y0DF1DX1.0 
Y0DF2  =  YGDF2DX1.0 
50  CONTINUE 
RETURN 

SUBROUTINE  IMPEDOC OMEGA, Z,ZR,ZI ,ZDB) 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

cx  THE  PURPOSE  OF  THIS  SUBROUTINE  IS  TO  CALCULATE  THE  IMPEDANCE, 

CX  AT  THE  DRIVING  POINT,  OF  A  CIRCULAR  DISK  BASED  UPON  BESSEL 

Cx  FUNCTIONS,  DISK  GEOMETRY,  AND  DRIVING  FREQUENCY. 

CX  xxxx  SPECIAL  CASE  OF  N=0  xxxx 
Cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


Cx 

LIST  OF 

VARIABLE 

X 

Cx 

ZZ 

COMPLEX  IMPEDANCE,  LBF-SEC/IN 

X 

Cx 

Z 

MAGNITUDE  OF  ZZ,  LBF-SEC/IN 

X 

Cx 

ZR 

REAL  PART  OF  ZZ,  LBF-SEC/IN 

X 

Cx 

ZI 

IMAGINARY  PART  OF  ZZ,  LBF-SEC/IN 

X 

Cx 

ZDB 

Z  IN  DB, 

X 

Cx 

A 

PLATE  RADIUS,  IN 

X 

Cx 

B 

WAVEGUIDE  ATTACHMENT  RADIUS,  IN 

X 

cx 

H 

PLATE  THICKNESS,  IN 

X 

Cx 

IMAT 

TYPE  OF  MATERIAL 

X 

Cx 

E 

YOUNG'S  MODULUS,  PSI 

X 

Cx 

G 

SHEAR  MODULUS,  PSI 

X 

CX 

D 

FLEXURAL  RIGIDITY,  PSI-IN 

X 

CX 

MU 

POISSON'S  RATIO 

X 

Cx 

RHO 

MASS  DENSITY  OF  PLATE,  SLUGS/INXX3 

X 

CX 

K 

SHEAR  COEFFICIENT 

X 

Cx 

OMEGA 

CIRCULAR  FREQUENCY  OF  INTEREST,  RAD/SEC 

X 

CX  ' 

DELTA1 

DEMENSIONL ESS  PARAMETER 

X 

Cx 

DELTA2 

DEMENSIONLESS  PARAMETER 

X 

CX 

DEL  TA3 

DEMENSIONLESS  PARAMETER 

X 

CX 

SIGMA1 

DEMENSIONLESS  PARAMETER 

X 

Cx 

SIGMA2 

DEMENSIONLESS  PARAMETER 

X 

CX 

RR 

DEMENSIONLESS  PARAMETER 

X 

Cx 

S 

DEMENSIONLESS  PARAMETER 

X 

CX 

L4 

DEMENSIONLESS  PARAMETER 

X 

CX 

PI 

3.14 - 

X 

CX 

ANSX 

INTERMEDIATE  STEP  X  ANSWERS 

X 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

IMPLICIT  C0MPLEXX8  (A,D,J,S,Y) 

C0MPLEXX8  ZZ,E,G 
REAL  MU , A , K 

COMMON/BESS1A/J1A, J1DF1A, J1DF2A, Y1A, Y1DF1A, Y1DF2A 


125 


*  *  tt  * 


C0MM0N/BESS1B/J1B, J1DF1B, J1DF2B, Y1B, Y1DF1B, Y1DF2B 
C0MM0N/BESS2A/J2A, J2DF1A, J2DF2A, Y2A, Y2DF1A, Y2DF2A 
COMMON/ BESS2B/ J 2B, J2DF1 B, J 2 DF2B, Y2B,Y2DF1B,Y2DF2B 
C0MM0N/BESS3A/J3A, J3DF1A, J3DF2A, Y3A, Y3DF1A, Y3DF2A 
C0MM0N/BESS3B/J3B, J3DF1B, J3DF2B, Y3B, Y3DF1B, Y3DF2B 
COMMON/GEOMTY/A , B , H , K 
COMMON/MAT/ I MAT , E, G, MU, RHO, TEMP 

COMMON/CONST1/DELTA1 , DELTA2, DELTA3, SIGMA1 , SIGMA2 
PI=ACOS( -1 . 0 ) 

IMPEDANCE  CALCULATION,  FROM  REDUCE  PROGRAM 

ANS4=Y2DF1BXJ2DF1A*Y1DF1A*J1DF1B*MUKSIGMA2*X2- 

2.*Y2DF1B*J2DF1AXY1DF1AXJ1DF1B*MU*SIGMA2X 

SIGMA1+Y2DF1BXJ2DF1AXY1DF1AXJ1DF1BXMUXSIGMA1 

XX2+Y2DF1BXJ2DF1AXYIDF1BXJ1DFIAXSIGMA2XX2X 

SIGMA1XDELTA1-Y2DF1BXJ2DF1AXY1DF1BXJ1DFIAX 

SIGMA2XX2XDELTA1-Y2DF1BXJ2DF1AXY1DF1BXJ1DF1AX 

SIGMA2XSIGMA1XX2XDELTA1+Y2DF1BXJ2DF1AXY1DF1BX 

J1DF1AXSIGMA2XSIGMA1XDELTA1-Y2DF1BXJ2DF1AXY1DF1B 

XJ1DF1AXMUXSIGMA2XX2+2.XY2DF1BXJ2DF1AXY1DF1B 

XJ1DF1AXMUXSIGMA2XSIGMA1-Y2DF1BXJ2DF1AXYIDF1B 

XJ1DF1AXMUXSIGMA1XX2 

AN53=Y2DF1AXJ2DF1BXYIDF1BXJ1DF1AXSIGMA2XSIGMA1 
XX2XDELTA1-Y2DF1AXJ2DF1BXY1DF1BXJ1DF1AXSIGMA2 
XSIGMA1XDELTA1+Y2DF1AXJ2DF1BXY1DF1BXJ1DFIAXMU 
XSIGMA2XX2-2.XY2DF1AXJ2DF1BXY1DFIBXJ1DF1AXMUX 
SIGMA2XSIGMA1+Y2DF1 Ax J2DF1BXY1DF1BXJ1DF1 AXMUX 
SIGMA1XX2+Y2DF1BXJ2DF2AXY1DF1AXJ1DF1BXSIGMA2 
XX2XSIGMA1XDELTA2-Y2DF1BXJ2DF2AXY1DF1AXJ1DF1B 
XSIGMA2XSIGMA1XX2XDELTA2-Y2DF1BXJ2DF2AXY1DF1AX 
J1DF1BXSIGMA2XSIGMA1XDELTA2+Y2DF1BXJ2DF2AXY1DF1A 
XJ1DF1BXSIGMA1XX2XDELTA2-Y2DF1BXJ2DF2AXY1DF1B 
XJ1DF1AXSIGMA2XX2XSIGMA1XDELTA2+Y2DF1BXJ2DF2A 
XY1DFIBXJ1DF1AXSIGMA2XSIGMA1XX2XDELTA2+Y2DF1BX 
J2DF2AXY1DF1BXJ1DF1AXSIGMA2XSIGMA1XDELTA2-Y2DF1B 
XJ2DF2AXY1DF1BXJIDF1AXSIGMAIXX2XDELTA2-Y2DF1B 
XJ2DF1AXY1DF2AXJ1DF1BXSIGMA2XX2XSIGMA1XDELTA1 
+Y2DF1BXJ2DF1AXY1DF2AXJ1DF1BXSIGMA2XX2XDELTA1 
+Y2DF1BXJ2DF1AXY1DF2AXJ1DF1BXSIGMA2XSIGMA1XX2 
XDELTA1-Y2DF1BXJ2DF1AXY1DF2AXJ1DF1BXSIGMA2X 
SIGMA1XDELTA1+ANSA 

AflS2  =  -Y2DF2Axj2DFl  BXY1DF1 AXJ1DF1  BXSIGMA2»*2X 
SIGMA 1 xDEl T A2+Y2DF2AX J2DF1BXY1DF1 Ax J1DFIB* 
SIGMA2XSIGMA1XX2XDELTA2+Y2DF2AXJ2DF1BXY1DF1AX 
J1DF1BXSIGMA2XSIGMA1XDELTA2-Y2DF2AXJ2DF1BXY1DF1A 
XJ1DF1BXSIGMA1XX2XDELTA2+Y2DF2AXJ2DF1BXYIDF1B 
XJ1DF1AXSIGMA2XX2XSIGMA1XDELTA2-Y2DF2AXJ2DF1B 
XY1DF1BXJ1DF1AXSIGMA2XSIGMA1XX2XDELTA2-Y2DF2AX 
J2DF1BXY1DF1BXJ1DF1AXSIGMA2XSIGMA1XDELTA2+Y2DF2A 
XJ2DF1BXY1DF1BXJ1DF1AXSIGMA1XX2XDELTA2+Y2DF1A 
XJ2DF1BXY1DF2AXJ1DF1BXSIGMA2XX2XSIGMA1XDELTAI 
-Y2DF1AXJ2DF1BXY1DF2AXJ1DF1BXSIGMA2XX2XDELTA1 
-Y2DF1AXJ2DF1BXY1DF2AXJ1DF1BXSIGMA2XSIGMA1XX2 
XDELTA1+Y2DF1AXJ2DF1BXY1DF2AXJ1DF1BXSIGMA2X 
SIGMA1XDELTA1-Y2DF1AXJ2DF1BXY1DF1AXJ1DF1BXMUX 
SIGMA2XX2+2.XY2DF1AXJ2DF1BXY1DF1AXJ1DF1BXMUX 
SIGMA2XSIGMA1-Y2DF1AXJ2DF1BXY1DF1AXJ1DF1BXMUX 
SIGMA1XX2-Y2DF1AXJ2DF1BXY1DF1BXJ1DF1AXSIGMA2 
XX2XSIGMA1XDELTA1+Y2DF1AXJ2DF1BXY1DF1BXJ1DF1A 
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.  XSIGMA2**2*DELTA1+ANS3 
ANS1=2 . XBXGXHXKXX2XPI*DELTA2XDELTA1XANS2 
ANS13=Y2BXJ2DF1BXY1DF1AXJ1DF1AXSIGMA2XSIGMA1XX 
.  2XDELTA1XX2-Y2BXJ2DF1BXY1DF1AXJ1DF1AXSIGMA2X 
.  SIGMA1XDELTA1XX2-Y2BXJ2DF1BXY1DF1AXJ1DF1AX 
.  SIGMA1XX2XDELTA1XX2+Y2BXJ2DF1BXY1DF1AXJ1DF1AX 

.  SIGMA1XDELTA1XK2 

ANS12=-2.XY2BXJ2DF1AXY1DF2AXJ1DF1BXSIGMA2X 
.  SIGMA1XDELTA1XX2+Y2BXJ2DF1AXY1DF2AXJ1DF1BX 
.  SIGMA2XDELTA1XX2-Y2BXJ2DF1AXY1DF1AXJ1DF1BXMUX 
.  SIGMA2XSIGMA1XDELTA1+Y2BXJ2DF1AXY1DF1AXJ1DF1B 
.  XMUXSIGMA2XDELTA1+Y2BXJ2DF1AXY1DF1AXJ1DF1BXMU 
.  XSIGMA1XX2XDELTA1-Y2BXJ2DF1AXY1DF1AXJ1DF1BXMU 
.  XSIGMA1XDELTA1-Y2BXJ2DF1AXY1DF1BXJ1DF1AX 
.  SIGMA2XSIGMA1XX2XDELTA1XX2+2.XY2BXJ2DF1AXY1DF1BX 
.  J1DF1AXSIGMA2XSIGMA1XDELTA1XX2-Y2BXJ2DF1AXY1DF1B 
.  XJ1DF1AXSIGMA2XDELTA1XX2+Y2BXJ2DF1AXY1DF1B 
.  XJ1DF1AXMUXSIGMA2XSIGMA1XDELTA1-Y2BXJ2DF1AX 
.  Y1DF1B*J1DF1AXMUXSIGMA2XDELTA1-Y2BXJ2DF1AX 
.  Y1DF1BXJ1DF1AXMUXSIGMA1XK2XDELTA1+Y2BXJ2DF1AX 
.  Y1DF1B*J1DF1AXMUXSIGMA1XDELTA1-Y2BXJ2DF1BX 
.  Y1DF2AXJ1DF1AXSIGMA2XSIGMA1XX2XDELTA1XX2+Y2BX 
.  J2DF1BXY1DF2AXJ1DF1AX3IGMA2XSIGMA1XDELTA1XX2+Y2B 
.  XJ2DF1BXY1DF2AXJ1DF1AXSIGMA1XX2XDELTA1XX2-Y2B 
.  XJ2DF1BXY1DF2AXJ1DF1AXSIGMA1XDELTA1X*2+ANS13 
ANS11=-Y2DF1BXJ2BXY1DF2AXJ1DF1AXSIGMA2XSIGMA1X 
.  DELTA1XX2-Y2DF1BXJ2BXY1DF2AXJ1DF1AXSIGMA1XX2X 
.  DELTA1XX2+YZDF1BXJ2BXY1DF2AXJ1DF1AXSIGMA1X 
.  DELTA1XX2-Y2DF1BXJ2BXY1DF1AXJ1DF1AXSIGMA2X 
.  SIGMA1XX2XDELTA1XX2+Y2DF1BXJ2BXY1DF1AXJ1DF1AX 
.  SIGMA2XSIGMA1XDELTA1XX2+Y2DF1BXJ2BXY1DF1AXJ1DF1A 
.  XSIGMA1XX2XDELTA1XX2-Y2DF1BXJ2BXY1DF1AXJ1DF1A 
.  XSIGMA1XDELTA1XX2-Y2BXJ2DF2AXY1DF1AXJ1DF1B 
.  XSIGMA2XSIGMA1XX2XDELTA2XDELTA1+Y2BXJ2DF2AXY1DF1A 
.  XJ1DF1BXSIGMA2XSIGMA1XDEITA2XDELTA1+Y2BXJ2DF2AX 
.  Y1DF1AXJ1DF1BXSIGMA1XX2XDELTA2XDELTA1-Y2BXJ2DF2A 
.  XY1DF1AXJ1DF1BXSIGMA1XDELTA2XDELTA1+Y2BX 
.  J2DF2AXY1DF1BXJ1DF1AXSIGMA2XSIGMA1XX2XDELTA2X 
.  DELTA1-Y2BX J2DF2AXY1DF1B* J1 DF1AXSIGMA2XSIGMA1 
.  XDELTA2XDELTA1-Y2BXJ2DF2AXY1DF1BXJ1DF1AX 
.  SIGMA1XX2XDELTA2XDELTA1+Y2BXJ2DF2AXY1DF1BXJ1DF1A 
.  XSIGMA1XDELTA2XDELTA1+Y2BXJ2DF1AXY1DF2AX 
.  J1DF1BXSIGMA2XSIGMA1XX2XDELTA1XX2+ANS12 
AM510=-Y2DF1BXJ2DF1AXY1DF2AXJ1BXSIGMA2XX2X 
.  SIGMA1XDELTA2XDELTA1+Y2DF1BXJ2DF1AXY1DF2AXJ1B 
.  XSIGMA2XX2XDELTA2XDELTA1+Y2DF1BXJ2DF1AXY1DF2AX 

.  J1BKSIGMA2XSIGMA1XDELTA2*DELTA1-Y2DF1BXJ2DF1AX 
.  Y1DF2AXJ1BXSIGMA2XDELTA2XDELTA1+Y2DF1BXJ2DF1A 
.  XY1DF1AXJ1BXMUXSIGMA2XX2XDELTA2-Y2DF1BXJ2DF1A 
.  XY1DF1AXJ1BXMUXSIGMA2XSIGMA1XDELTA2-Y2DF1BX 
.  J2DF1AXY1DF1AXJ1BXMUXSIGMA2XDELTA2+Y2DF1BX 
.  J2DF1AXY1DF1AXJ1BXMUXSIGMA1XDELTA2+Y2DF1BX 
.  J2DF1AXY1BXJ1DF1AXSIGMA2XX2XSIGMA1XDELTA2XDELTA1 
.  -Y2DF1BXJ2DF1AXY1BXJ1DF1AXSIGMA2XX2XDELTA2X 
.  DELTA1-Y2DF1BXJ2DF1AXY1BXJ1DF1AXSIGMA2XSIGMA1 
.  XDELTA2XDELTA1+Y2DF1BXJ2DF1AXY1BXJ1DFIAX 
.  SIGMA2XDELTA2XDELTA1-Y2DF1BXJ2DF1AXY1BXJ1DF1A 
.  XMUX5IGMA2XX2XDELTA2+Y2DF1BXJ2DF1AXY1BXJ1DF1A 
.  XMUXSIGMA2XSIGMA1XDELTA2+Y2DF1BXJ2DF1AXY1BX 
.  J1DF1AXMUXSIGMA2XDELTA2-Y2DF1BXJ2DF1AXY1BX 
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J1DF1A*MU*SIGMA1*DELTA2+Y2DF1B*J2B*Y1DF2AX 

J1DF1Ax$IGMA2xSIGMA1**2*DELTA1**2+ANS11 

ANS9=Y2DF1AXJ2B*Y1DF1AXJ1DF1B*MU*SIGMA2XSIGMA1 
*DELTA1-Y2DF1AXJ2B*Y1DF1AXJ1DF1B*MUXSIGMA2X 
DELTA1-Y2DF1AXJ2BXY1DF1AXJ1DF1BXMUXSIGMA1XX2X 
DELTA1+Y2DF1AXJ2BXY1DF1AXJ1DF1BXMUXSIGMA1X 
DELTA1+Y2DF1AXJ2BXY1DF1BXJ1DF1AXSIGMA2XSIGMA1 
XX2XDELTA1XX2-2.XY2DF1AXJ2BXY1DF1BXJ1DF1AX 
SIGMA2XSIGMA1XDELTA1XX2+Y2DF1AXJ2BXY1DF1BXJ1DF1A 
XSIGMA2XDELTA1XX2-Y2DF1AXJ2BXY1DF1BXJ1DF1A 
XMUXSIGMA2XSIGMA1XDELTA1+Y2DF1AXJ2BXY1DF1BX 
J1DF1AXMUXSIGMA2XDELTA1+Y2DF1AXJ2BXY1DF1BX 
J1DF1AXMUXSIGMA1XX2XDELTA1-Y2DF1AXJ2BXY1DF1BX 
J1DF1AXMUXSIGMA1XDELTA1+Y2DF1BXJ2DF2AXY1DF1AX 
J1BXSIGMA2XX2XSIGMA1XDELTA2XX2-2 . XY2DF1BXJ2DF2AX 
Y1DF1AXJ1BXSIGMA2XSIGMA1XDELTA2XX2+Y2DF1B*J2DF2A 
XY1DF1AXJ1BXSIGMA1XDELTA2XX2-Y2DF1BXJ2DF2A 
XY1BXJ1DF1AXSIGMA2XX2XSIGMA1XDELTA2XX2+2 . XY2DF1B 
XJ2DF2AXY1BXJ1DF1AXSIGMA2XSIGMA1XDELTA2XX2- 
Y2DF1BXJ2DF2AXY1BXJ1DF1AXSIGMA1XDELTA2XX2+ 

ANSI  0 

ANS8*-Y2DF1AXJ2DF1BXY1DF2AXJ1BXSIGMA2XSIGMA1X 

DELTA2XDELTA1+Y2DF1AXJ2DF1BXY1DF2AXJ1BXSIGMA2 

XDELTA2XDELTA1-Y2DF1AXJ2DF1BXY1DF1AXJ1BXMUX 

SIGMA2XX2XDELTA2+Y2DF1AXJ2DF1BXY1DF1AXJ1BXMUX 

SIGMA2XSIGMA1XDELTA2+Y2DF1AXJ2DF1BXY1DF1AXJ1B 

XMUXSIGMA2XDELTA2-Y2DF1AXJ2DF1BXY1DF1AXJ1BXMU 

XSIGMA1XDELTA2-Y2DF1AXJ2DF1BXY1BXJ1DF1AX 

SIGMA2XX2XSIGMA1XDELTA2XDELTA1+Y2DF1AXJ2DF1BXY1B 

XJ1DF1AXSIGMA2XX2XDELTA2XDELTA1+Y2DF1AXJ2DF1BX 

Y1BXJ1DF1AXSIGMA2XSIGMA1XDELTA2XDELTA1-Y2DF1AX 

J2DF1BXY1BXJ1DF1AXSIGMA2XDELTA2XDELTA1+Y2DF1A 

XJ2DF1BXY1BXJ1DF1AXMUXSIGMA2XX2XDELTA2-Y2DF1A 

XJ2DF1BXY1BXJ1DF1AXMUXSIGMA2XSIGMA1XDELTA2- 

Y2DF1AXJ2DF1BXY1BXJ1DF1AXMUXSIGMA2XDELTA2+ 

Y2DF1AXJ2DF1BXY1BXJ1DF1AXMUXSIGMA1XDELTA2- 

Y2DF1AXJ2BXY1DF2AXJ1DF1BXSIGMA2XSIGMA1XX2X 

DELTA1XX2+2.XY2DF1AXJ2BXY1DF2AXJ1DF1BXSIGMA2X 

SIGMA1XDELTA1XX2-Y2DF1AXJ2BXY1DF2AXJ1DF1BX 

SIGMA2XDELTA1XX2+ANS9 

AHS7=-Y2DF2AXJ2BXY1DF1AXJ1DF1BXSIGMA1XX2X 
DELTA2*DELTA1+Y2DF2AXJ2B*Y1DF1 AXJ1DF1 BXSIGMA1 
XDELTA2XDELTA1-Y2DF2AXJ2BXY1DF1BXJ1DF1AX 
SIGMA2XSIGMA1XX2XDELTA2XDELTA1+Y2DF2AXJ2BXY1DF1B 
XJ1DF1AXSIGMA2XSIGMA1XDELTA2XDELTA1+Y2DF2AXJ2BX 
Y1DF1BXJ1DF1AXSIGMA1XX2XDELTA2XDELTA1-Y2DF2AXJ2B 
XY1DF1BXJ1DF1AXSIGMA1XDELTA2X0ELTA1-Y2DF1AX 
J2DF2AXY1DF1BXJ1BXSIGMA2XX2XSIGMA1XDELTA2XX2+ 
Y2DF1AXJ2DF2AXY1DF1BXJ1BXSIGMA2XX2XDELTA2XX2+ 
Y2DF1AXJ2DF2AXY1DF1BXJ1BXSIG!'JA2XSIGMA1XDELTA2 
XX2-Y2DF1AXJ2DF2AXY1DF1BXJ1BXSIGMA2XDELTA2XX2 
+Y2DF1AXJ2DF2AXY1BXJ1DF1BXSIGMA2XX2XS1GMA1X 
DELTA2XX2-Y2DF1AXJ2DF2AXY1BXJ1DF1BXSIGMA2XX2X 
DELTA2XX2-Y2DF1AXJ2DF2AXY1BXJ1DF1BXSIGMA2X 
SIGMA1XDELTA2XX2+Y2DF1AXJ2DF2AXY1BXJ1DF1BX 
SIGMA2XDELTA2XX2+Y2DF1AXJ2DF1BXY1DF2AXJ1BX 
SIGMA2XX2XSIGMA1XDELTA2XDELTA1-Y2DF1AXJ2DF1BXY1DF2A 
XJ1BXSIGMA2XX2XDELTA2XDELTA1+ANS8 

ANS6=Y2DF2Axj2DF1A»Y1DF1BxJ1BxSIGMA2xx2xSIGMA1 

XDEITA2XX2-Y2DF2AXJ2DF1AXY1DF1BXJ1BXSIGMA2XX2 
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non  noo  non  non  non 


.  XDELTA2**2-Y2DF2A*J2DFlA*yiDFlBXJlBX$IGMA2* 

.  SIGMA1XDELTA2XX2+Y2DF2AXJ2DF1AXY1DF1BXJ1BX 
.  SIGMA2XDELTA2XX2-Y2DF2AXJ2DF1AXY1BXJ1DF1BX 
.  SIGMA2XX2XSIGMA1XDELTA2XX2+Y2DF2AXJ2DF1AXY1BX 
.  J1DF1BXSIGMA2XX2XDELTA2XX2+Y2DF2AXJ2DF1AXY1BX 
.  J1DF1BXSIGMA2XSIGMA1XDELTA2XX2-Y2DF2AXJ2DF1AXY1B 
.  XJ1DF1BXSIGMA2XDELTA2XX2-Y2DF2AXJ2DF1BXY1DF1A 
.  XJ1BXSIGMA2XX2XSIGMA1XDELTA2XX2+2.XY2DF2AXJ2DF1B 
.  XY1DF1AXJ1BXSIGMA2XSIGMA1XDELTA2XX2-Y2DF2AX 
.  J2DF1BXY1DF1AXJ1BXSIGMA1XDELTA2XX2+Y2DF2AX 
.  J2DF1BXY1BXJ1DF1AXSIGMA2XX2XSIGMA1XDELTA2XX2-2.X 

.  Y2DF2AXJ2DF1BXY1BXJ1DF1AXSIGMA2XSIGMA1XDELTA2 
.  XX2+Y2DF2AXJ2DF1BXY1BXJ1DF1AXSIGMA1XDELTA2XX2 
.  +Y2DF2AXJ2BXY1DF1AXJ1DF1BXSIGMA2XSIGMA1XX2X 
.  DELTA2XDELTA1-Y2DF2AXJ2BXY1DF1AXJ1DF1BXSIGMA2 
.  XSIGMA1XDELTA2XDELTA1+ANS7 
ANS5=AX0MEGAXANS6X( 0 . 0 , -1 . 0) 

IMPEDANCE 

ZZ  =  ANS1/ANS5 

REAL  PART  OF  IMPEDANCE 

ZR  =  REAL(ZZ) 

IMAGINARY  PART  OF  IMPEDANCE 

ZI  =  AIMAG(ZZ) 

MAGNITUDE  OF  IMPEDANCE 

Z  =  CABS(ZZ) 

DB  OF  IMPEDANCE 

ZDB  =  20 . OXLOGIO(Z) 

RETURN 

END 
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Figure  66.  Experimental  Real  Part  of  the  Driving  Point  Impedance  of 
a  6  in  Radius  Elastic  Hate,  With  and  Without  Shaker  and 
Mount  Correction,  in  the  Frequency  Range  of  5  to  106  Hz. 


Figure  67.  Experimental  Imaginary  Part  of  the  Driving  Point 
impoHanfo  of  a  6  in  Radius  Elastic  Hate,  With  and 
Without  fihakflr  and  Mount  Correction,  in  the  Frequency 
Range  of  5  to  105  Hz. 
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Figure  68.  Experimental  Real  Part  of  the  Driving  Point  Impedance  of 
a  6  in  Radius  Viscoelastic  Plate,  With  and  Without  Shafaw 
and  Mount  Correction,  in  the  Frequency  Range  of  100  to 
2000  Hz  at  a  Temperature  of  100J5  Deg.  F. 
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Figure  69.  Experimental  Imaginary  Part  of  the  Driving  Point 

impedance  of  a  6  in  Radius  Viscoelastic  Plate,  With  and 
Without  Shaker  and  Mount  Connection,  in  the  Frequency 
Range  of  100  to  2000  Hz  at  a  Temperature  of  100l5 
Deg.  F. 
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Figure  70.  Effect  of  Shear  Coefficient  on  the  Theoretical  Real  Part  of 
the  Driving  Point  Impedance  for  a  5  in  Radius  Elastic 
Plate  in  the  Frequency  Range  of  10  to  106  Hz. 
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Figure  7L  Effect  of  Shear  Coefficient  on  the  Theoretical  Imaginary 
Part  of  the  Driving  Point  Impedance  for  a  5  in  Radius 
Elastic  Plate  in  the  Frequency  Range  of  10  to  105  Hz. 
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Figure  72. 


Effect  of  Shear  Coefficient  on  the  Theoretical  Real  Part  of 
the  Driving  Point  Impedance  for  a  5  in  Radius  Elastic 
Plate  in  the  Frequency  Range  of  100  to  2000  Hz. 
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Figure  73.  Effect  of  Shear  Coefficient  cm  the  Theoretical  Imaginary 
Part  of  the  Driving  Point  Impedance  for  a  5  in  Radius 
Elastic  Plate  in  the  Frequency  Range  of  100  to  2000  Hz. 
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Figure  75. 
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Effect  of  Shear  Coefficient  and  Material  Characteristics  on 
the  Theoretical  Imaginary  Part  of  the  Driving  Point 
Impedance  for  a  5  in  Radius  Viscoelastic  Plate  in  the 
Frequency  Range  of  10  to  105  Hz  at  a  Temperature  of  80.0 
Deg.  F. 
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Figure  77.  Effect  of  Shear  Coefficient  and  Material  Characteristics  on 
the  Theoretical  Imaginary  Part  of  the  Driving  Point 
Impedance  for  a  5  in  Radius  Viscoelastic  Plate  in  the 
Frequency  Range  of  100  to  2000  Hz  at  a  Temperature  of 
77.3  Deg.  F. 
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Figure  78.  Experimental  Real  Part  of  Hie  Driving  Point  Impedance  of 
a  5  in  Radius  Elastic  Plate,  With  and  Without  Shaker  and 
Mount  Correction,  in  the  Frequency  Range  of  5  to  106  Hz. 
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Figure  79.  Experimental  Imaginary  Part  of  the  Driving  Point 
Impedance  of  a  5  in  Radius  Elastic  Plate,  With  and 
Without  Shaker  and  Mount  Correction,  in  the  Frequency 
Range  of  5  to  105  Hz. 


80.  Experimental  Real  Part  of  the  Driving  Point  Impedance  of 
a  5  in  Radius  Viscoelastic  Plate,  With  and  Without  Shaker 
and  Mount  Correction,  in  the  Frequency  Range  of  5  to  105 
Hz  at  a  Temperature  of  80.0  Deg.  F. 
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Figure  8L  Experimental  Imaginary  Part  at  Hie  Driving  Point  f 

Impedance  of  a  5  in  Radius  Viscoelastic  Plate,  With  and 
Without  Shaker  and  Mount  Conection,  in  the  Frequency 
Range  of  5  to  105  Hz  at  a  Temperature  of  80.0  Deg.  F. 

I 

146 


I 


oos 


ai 

Ed 

Gg 

CiCd 

oe 

zzo 

K30 

g-fJ! 

_3  <r<; 

zz 

uu 


CEOi 

CdCd 

O.O. 

XX 

ww; 

— rvi 


O’Ot  O'OG  002  001 

(NI/03S-JS1)  33MVCI3dKI 


Figure  82.  Experimental  Real  Part  of  the  Driving  Point  Impedance  of 
a  5  in  Radius  Viscoelastic  Plate,  With  and  Without  Shaker 
and  Mount  Correction,  in  the  Frequency  Range  of  100  to 
2000  Hz  at  a  Temperature  of  77.3  Deg.  F. 


146 


Figure  83.  Experimental  Imaginary  Part  of  die  Driving  Paint 

Impedance  of  a  5  in  Radius  Viscoelastic  Plate,  With  anH 
Without  Shaker  and  Mount  Correction,  in  die  Frequency 
Range  of  100  to  2000  Hz  at  a  Temperature  of  77.3  Deg.  F. 
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Figure  85,  Theoretical  and  Experimental  Comparison  of  the 

Imaginary  Part  of  the  Driving  Point  Impedance  for  a  5  in 
Radius  Elastic  Hate  in  the  Frequency  Rang*  of  5  to 
105  Hz. 
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Figure  86.  Theoretical  and  Experimental  Comparison  of  the  Real  Part 
of  the  Driving  Point  Impedance  for  a  6  in  Radius  Elastic 
Plate  in  the  Frequency  Range  of  100  to  2000  Hz. 
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Figure  87.  Theoretical  and  Experimental  Comparison  of  the 

Imaginary  Part  of  the  Driving  Point  Impedance  for  a  5  in 
Radius  Elastic  Hate  in  the  Frequency  Range  of  100  to 
2000  Hz. 
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Figure  88.  Theoretical  and  Experimental  Comparison  of  the  Real  Part 
of  the  Driving  Point  Impedance  for  a  5  in  Radius 
Viscoelastic  Hate  in  the  Frequency  Range  of  5  to  105  Hz  at 
a  Temperature  of  80.0  Deg.  F. 
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Figure  89.  Theoretical  and  Experimental  Comparison  of  the 


Imaginary  Part  of  the  Driving  Point  Impedance  for  a  5  in 
Radius  Viscoelastic  Plate  in  the  Frequency  Range  of  5  to 
105  Hz  at  a  Temperature  of  80.0  Deg  F. 
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l.  Theoretical  and  Experimental  Comparison  of  the  Real  Part 
of  the  Driving  Point  Impedance  for  a  5  in  Radius 
Viscoelastic  Plate  in  the  Frequency  Range  of  100  to  2000  Hz 
at  a  Temperature  of  77.3  Deg.  F. 
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Figure  9L  Theoretical  and  Experimental  Comparison  of  the 

Imaginary  Part  of  the  Driving  Point  Impedance  for  a  5  in 
Radius  Viscoelastic  Plate  in  die  Frequency  Range  of  100  to 
2000  Hz  at  a  Temperature  of  77.3  Deg.  F. 
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Figure  92.  Effect  of  Shear  Coefficient  on  the  Theoretical  Real  Part  of 
the  Driving  Point  Impedance  for  a  4  in  Radius 
Rate  in  the  Frequency  Range  of  10  to  106  Hz. 
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Figure  93.  Effect  of  Shear  Coefficient  on  the  Theoretical  Imaginary 
Part  of  the  Driving  Point  Impedance  for  a  4  in  Radius 
FTInstir  Plate  in  the  Frequency  Range  of  10  to  106  Hz. 
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Figure  94.  Effect  of  Shear  Coefficient  on  the  Theoretical  Real  Part  of 
the  Driving  Point  Impedance  for  a  4  in  Radius  Elastic 
Plate  in  the  Frequency  Range  of  100  to  2000  Hz. 
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Figure  95.  Effect  of  Shear  Coefficient  on  the  Theoretical  Imaginary 
Part  of  the  Driving  Point  Impedance  for  a  4  in  Radius 
Elastic  Hate  in  the  Frequency  Range  of  100  to  2000  Hz. 


.A  r 


Figure  96. 
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Effect  of  Sh^nr  Coefficient  and  Material  Characteristics  on 
the  Theoretical  Real  Part  of  the  Driving  Point  Impedance 
for  a  4  in  Radius  Viscoelastic  Plate  in  the  Frequency 
Range  of  10  to  105  Hz  at  a  Temperature  of  7&5  Deg.  F. 
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Figure  97.  Effect  of  Shear  Coefficient  and  Material  Characteristics  on 
the  Theoretical  Imaginary  Part  of  the  Driving  Point 
Impedance  for  a  4  in  Radius  Viscoelastic  Plate  in  the 
Frequency  Range  of  10  to  105  Hz  at  a  Temperature  of 
78.5  Deg.  F. 
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Figure  10L  Experimental  Imaginary  Part  of  the  Driving  Point 
Impedance  of  a  4  in  Radius  Elastic  Plate,  With  and 
Without  Shaker  and  Mount  Correction,  in  the  Frequency 
Range  of  5  to  105  Hz. 
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Figure  102.  Experimental  Real  Part  of  the  Driving  Point  Impedance  of 
a  4  in  Radius  Viscoelastic  Plate,  With  and  Without  Shaker 
and  Mount  Correction,  in  the  Frequency  Range  of  5  to  105 
Hz  at  a  Temperature  of  7&5  Deg.  F. 
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Figure  103.  Experimental  Imaginary  Part  of  the  Driving  Point 

Impedance  of  a  4  in  Radius  Viscoelastic  Plate,  With  nnH 
Without  Shaker  and  Mount  Correction,  in  the  Frequency 
Range  of  5  to  105  Hz  at  a  Temperature  of  78J5  Deg.  F. 
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Figure  104^  Experimental  Real  Part  o f  ftw  Driving  Point  Impedance  of 
a  4  in  Radius  Viscoelastic  Plate,  With  Without  Shaker 
and  Mount  Correction,  in  the  Frequency  Range  of  100  to 
2000  Hz  at  a  Temperature  of  74J5  Deg.  F. 
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Figure  105.  Experimental  Imaginary  Part  of  the  Driving  Point 

Impedance  of  a  4  in  Radios  Viscoelastic  Plate,  With  and 
Without  SHflkftr  and  Mount  Correction,  in  the  Frequency 
Range  of  100  to  2000  Hz  at  a  Temperature  of  74.5  Deg.  F. 
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Figure  106.  Theoretical  and  Experimental  Comparison  of  the  Real  Part 
of  the  Driving  Point  Impedance  foo*  a  4  in  Radius  Elastic 
Plate  in  the  Frequency  Range  of  5  to  105  Hz. 
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Figure  107.  Theoretical  and  Experimental  Comparison  of  the 


Imaginary  Part  of  the  Driving  Point  Impedance  for  a  4  in 
Radius  Elastic  Plate  in  the  Frequency  Range  of  5  to 
105  Hz. 
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Figure  108,  Theoretical  and  Experimental  Comparison  of  the  Real  Part 
of  the  Driving  Point  Impedance  for  a  4  in  Radius  Elastic 
Plate  in  the  Frequency  Range  of  100  to  2000  Hz. 
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Figure  109.  Theoretical  and  Experimental  Comparison  of  the 

Imaginary  Part  of  the  Driving  Point  Impedance  for  a  4  in 
Radius  Elastic  Plate  in  the  Frequency  Range  of  100  to 


Figure  110.  Theoretical  and  Experimental  Comparison  of  the  Real  Part 
of  the  Driving  Point  Impedance  for  a  4  in  Radius 
Viscoelastic  Plate  in  the  Frequency  Range  of  5  to  106  Hz  at 


a  Temperature  of  7&5  Deg.  F. 
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Figure  11L  Theoretical  and  Experimental  Comparison  of  the 

Imaginary  Part  of  the  Driving  Point  Impedance  for  a  4  in 
Radius  Viscoelastic  Plate  in  the  Frequency  Range  of  5  to 
105  Hz  at  a  Temperature  of  7&5  Deg.  F. 
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Figure  112-  Theoretical  and  Experimental  Comparison  of  the  Real  Part 
of  the  Driving  Point  Impedance  for  a  4  in  Radius 
Viscoelastic  Hate  in  the  Frequency  Range  of  100  to  2000  Ha 
at  a  Temperature  of  74J5  Deg.  F. 
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